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ABSTRACT 


‘A  method  for  calculating  the  response  of  simply  supported  windows  to 
sonic  booms  haa  been  developed.  The  procedure  is  based  on  a  linear  one- 
degree-of-freedom  analysis  plus  estimates  of  the  importance  of  nonlinear 
and  multimodal  effects.  Effects  of  stress  raisers  and  of  movement  followed 
by  impact  of  loose  windows  are  not  considered. 

Significant  contributions  to  the  maximum  stress  in  windows  subjected 
to  2  psf  aonic  booms  are  made  by  large  deflections  (nonlinearities),  modes 
above  the  fundamental,  and  the  internal  pressure  built  up  in  the  building 
by  the  boom. 

An  attempt  to  estimate  statistically  the  occurrence  of  window  failure 
under  2  psf  booms  was  frustrated  by  the  lack  of  precise  knowledge  of  the 
statistical  distribution  of  glass  strength. 


Hi 


PREFACE 


This  report  ia  one  of  a  aeriea  of  technical  reporta  dealing  with 
varioua  effects  of  aonic  booms.  The  research  waa  sponsored  by  tho 
National  Sonic  Boom  Evaluation  Office  of  the  Air  Force. 
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RESPONSE  OF  WINDOWS  TO  SONIC  BOOHS 


SECTION  I 
INTRODUCTION 

The  reaponse  of  windows  to  aonic  hoorn  loadings  preaenta  an  impor¬ 
tant  prohlem  in  the  evaluation  of  sonic  hoorn  effects  on  structures.  Ex¬ 
perimental  studios  which  have  heen  undertaken  to  ascertain  the  effects 
of  hooma  on  windows  of  various  sizes  are  those  of  Maglieri,  Huckel  and 
Perrott  (1961),  Blume  and  Associates  (1965),  and  Freynik  (1963).  These 
teat  progrems  have  ahown  that  nearly  all  windows  ere  hroken  hy  hoorn 
preaaurea  of  100  paf  and  that  e  few  will  hreek  at  preaaures  of  10  or  20 
psf.  No  lower  limit  has  heen  estahlished  on  the  pressure  required  to 
hreek  a  window.  Theae  three  studies  and  the  stetic  teata  of  Bowles  end 
Sugsrman  (1952)  and  Orr  (1957)  heve  shown  that  windows  have  a  nonlinear 
hehavior  even  at  low  preaaurea. 

Nonlinearities  hecotne  important  when  the  central  deflection  exceeds 
one-half  the  plate  thickneaa,  and,  for  uaual  window  dimensions,  static 
preaaurea  of  1  paf  to  10  paf  will  cause  auch  a  deflection.  The  conse¬ 
quence  of  the  nonlinearitiea  ia  that  atreaa  and  deflection  ere  not  pro¬ 
portional  to  pressure  or  to  each  other.  Hence,  there  has  heen  consider¬ 
able  difficulty  in  interpreting  experimental  reaulta  and  in  extrapolating 
the  findings.  In  moat  esaes  the  data  were  not  correlated  with  theory. 

In  the  dynamic  teata  the  pressures  acting  on  the  window  were  not  meeaured, 
only  the  nominal  preaaure  in  the  vicinity  waa  measured.  In  esaes  where 
window  deflection  was  measured,  only  maximum  excursion  waa  ohtained. 

Knowledge  of  window  hehavior  requires  careful  procuring  of  the  right 
data  and  data  reduction  hesed  on  3  theoretical  analyais  of  the  window 
motion.  This  present  study  ia  intended  to  hring  together  available  theo¬ 
retical  knowledge  on  the  suhject  to 

1.  indicete  the  experimental  parameters  which  ahould  he  meesured, 

2.  provide  a  baais  for  data  reduction, 

3.  provide  a  haais  for  prediction  of  window  atreasea  and  deflec¬ 
tions  in  response  to  booms, 
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4.  indicata  those  araas  whara  furthar  axparimental  or  theoratlcal 
work  Is  moat  naadad. 

Tho  study  is  primarily  concerned  with  boom  loadings  with  peak  presaures 
of  about  2  psf . 

Nat  loading  on  a  window  is  a  function  of  tha  stiffneas  and  volume 
of  tha  structure  in  which  tha  window  occurs.  Ra^ponse  of  a  window  to 
loading  is  highly  nonlinear  and  dependa  on  the  participation  of  savaral 
deflection  modes.  Because  of  this  complexity  tha  problem  has  not  baen 
aolvad  analytically  but  has  been  approximated  by  considering  (l)  general 
properties  of  windowa,  (2)  the  linear  responsa  of  windows  in  the  funda¬ 
mental  mode,  (3)  linear  responsa  in  all  modes,  (4)  effect  of  internal 
preasure  on  the  reaponae,  and  (5)  nonlinear  reaponsa  in  the  first  moda  . 
By  comparing  tha  results  of  (2),  (3),  (4),  and  (5),  wa  have  estimated 
the  nonlinear,  multimodal  reaponsa.  Finally,  using  tha  calculated  res- 
ponae,  the  atatiatical  probability  of  window  damaga  was  estimated. 


SECTION  II 


GENERAL  NATURE  OF  WINDOWS 

To  celculste  window  response,  we  need  to  know  the  normal  mounting 
conditions  of  windows,  strength  snd  modulus  of  glass,  common  sires  of 
windows,  end  the  neturel  iroquoncies  of  windows. 

Window  glsss  is  mounted  either  with  glazing  nsils  snd  putty  or  in 
s  rubber  sesl.  It  may  be  assumed  that  the  edge  conditions  cen  be  approxi¬ 
mated  by  s  simple  support  with  e  moveable  edge.  This  sssumpticn  hss  been 
studied  by  Freynlk  (1963).  Freynik  used  square  windows  mounted  with 
glezing  nsils  snd  putty.  He  compared  empirical  nsturel  frequencies  with 
those  from  violation  theory  for  simple  supported  pistes,  snd  compared 
measured  stresses  with  theoretical  vslues.  Goth  comparisons  showed  thst 
the  experimental  results  could  be  explained  on  the  bssls  of  s  simple 
support  condition.  The  results  of  Bowles  snd  Sugsrmsn  (1952)  in  stetlc 
testing  of  windows  confirms  the  existence  of  the  simple,  moveable  boun¬ 
dary  condition. 

An  unknown,  but  probably  large,  fraction  of  existing  windows  sre 
loose  in  their  mountings  snd/or  have  stress  veisers  present  in  poor  mount¬ 
ings.  The  results  of  this  report  csnnot  be  applied  to  such  defective 
windows  since  the  sbove  assumption  of  simple  support  edge  conditions  does 
not  spply.  Our  results  should  spply  well  to  calculation  of  the  response 
of  new,  well-mounted  windows  such  ss  those  usually  encountered  in  labora¬ 
tory  tests  or  field  test  structures. 

Strength  of  window  glsss  has  been  reported  by  R.  W.  McKinley  (1964); 
vslues  for  strength  snd  other  parameters  of  glass  windows  ere  lifted  in 
Table  1.  Strengths  were  determined  by  the  stsnderd  A STM  (1965)  beam  test 
procedure.  Flexursl  strength  was  evslueted  as  failure  in  tension:  glsss 
fells  without  sny  significant  amount  of  plsstic  deformation  so  thst  the 
yield  value  equals  the  ultimate  strength.  There  is  considersble  scatter 
in  results  from  strength  tests.  McKinley  (1964)  states  thst  the  strengths 
are  normally  distributed  snd  suggests  using  s  coefficient  of  variation  of 
25$  (that  is,  the  standard  deviation  is  one-fourth  the  mean  strength). 
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TABLE  I 


PROPERTIES  OF  WINDOW  GLASS 
(from  McKinley ) 


Property 

Velue 

Flexurel  Strength 

Short  duration:  sonic  booms,  blests 

One  minute  lo*ding:  wind 

Elestic  modulue 

Density 

Poisson’s  ratio 

6600  psi 

4400  psi 

107  ps  i 

0.09  lb/in3 

0.23 

On  the  other  hand;  Shand  (1958)  states  that  the  distribution  it  not  nor¬ 
mal,  but  skewed  so  thet  the  mean  is  larger  than  median  and  mode,-  possibly 
a  lognormel  or  Poisson  distribution  would  fit  the  data. 

Normal  sizes  for  glase  panes  are  dictated  by  building  code  require¬ 
ments  which  in  turn  are  based  on  wind  loadings  end  nominal  factors  of 
safety.  The  allowed  thicknesees  for  panes  of  various  ereas  as  specified 
in  the  Uniform  Building  Code  (1964)  are  shown  in  the  first  two  columns 
of  Teble  2.  The  values  are  for  a  wind  loading  of  20  psf.  Dimension  of 
the  window  may  be  combined  to  form  the  dimensionless  parameter,  a/h, 
where  a  is  the  length  of  one  side  of  a  square  pane  and  h  is  the  thickness , 
For  a  rectangular  pane,  e  is  taken  ss  the  square  root  of  the  area.  In  a 
equere  pane  the  a/h  ratio  and  the  pressure  govern  the  magnitude  of  stress 
in  the  pane  under  uniform  loading.  From  the  building  code  requirements, 
the  maximum  allowable  retlo  (a/h)  was  determined  and  listed  in  column  3 
of  Table  2.  If  it  is  presumed  that  a  builder  will  alwsys  use  the  minimum 
thickness  permitted,  then  we  cen  also  obtain  a  minimum  ratio  ss  in  column 
4  of  Table  2.  This  minimum  value  of  a/h  was  computed  for  each  thickness 
using  the  value  of  e  associated  with  the  next  smaller  thickness.  The 
minimum  values  are  listed  opposite  the  thickness  values  used  in  each  com¬ 
putation. 


TABLE  2 

Nominal  Sizes  of  Square  Panes* 


Thickness 

t, 

_  (in.) 

Maximum  Allowed 
a/h 

(dimensionless) 

Minimum  Probable 

e/h 

(dimensionless) 

5.8 

0.085 

340 

141** 

10.85 

0.115 

343 

251 

12 

i/e 

332 

316 

27 

3/ 16 

332 

222 

48 

7/32,  1/4 

380,  333 

285,  242 

75 

5/ 16 

333 

26B 

108 

3/8 

333 

277 

190 

1/2 

331 

250 

*  Data  in  Columns  (l)  end  (2)  ere  from  the  Uniform  Building  Code  (1964). 
**  Based  on  e  12  x  12-inch  window. 


The  listings  in  Table  2  show  that  the  range  of  a/h  values  ia  not 
very  large.  Since  most  of  the  values  are  between  220  and  340,  this  range 
wes  used  in  the  enelyaea  of  this  report. 

The  netural  frequencies  of  windows  ere  importent  for  studies  of  dy¬ 
namic  loadings.  For  this  analysis  the  windows  are  treated  as  simply 
supported  plates  undergoing  small  deflections.  This  epproach  appeers  to 
be  adequate  for  deflections  which  ere  leaa  than  the  pane  thickneua. 

Figure  1  ia  a  greph  of  natural  frequency  as  a  function  of  the  aree-to- 
thickneos  ratio,  calculated  from  the  equation  shown  on  the  figure.  Square 
end  rectengular  panes  with  severel  eapect  ratios  ere  included  end  the 
firat  three  frequencies  ere  considered.  These  three  frequencies  corres¬ 
pond  to  modes  in  which  there  are  one  maximum  point  of  deflection,  three 
maxima,  end  nine  maxime.  The  renge  of  erea-to-thickness  ratios  for  var¬ 
ious  glass  thicknesses  is  also  shown  to  indicate  the  probeble  frequencies 
which  are  encountered. 


5 


(ft1) 


krHMu 

.m 


it:/1!*  i 
:ksi.  —  ~ 


iniliitdyi 


",  •  .  nr”  1 

■«•«>  •  r 

"  III  ":i.  Ill  lilHU 

!l.''‘S!,ai':86!tiS:v 


•  <tVUlHNN<  f* 
IfRHDHElii  • 


... 

IJJf 
-  -  1  » 

'  !!H  . 

SECTION  III 


DEFLECTION  AND  STRESS  OF  WINDOWS  UNDER  SONIC  BOOHS 

A,  Problem  Statement 

To  determine  the  deflection  end  stress  of  windows  under  booms  let 
us  examine  quslitstively  pressures  on  the  buildings  end  reaction  of  the 
building.  Cheracteriatic  boom  pressures  related  to  structure  reection 
are  shown  in  Fig.  2.  In  the  free  field  tbe  pressure  rises  sharply.  At 
the  front  wall  of  the  building  the  pressure  shows  a  sharp  rise  and  evi¬ 
dence  of  reflections  from  tbe  ground  end  wall.  The  boom  pressure  at  the 
rear  shows  e  more  gradual  rise  and  folios  that  on  the  front  face  by  s 
few  milliseconds.  The  pressure  at  each  point  on  the  roof  has  s  sharp 
rise,  but  the  averege  roof  pressure  hta  a  longer  rise  time,  corresponding 
to  the  travel  time  of  the  wave  across  the  roof. 

The  effects  of  these  loadings  on  tbe  building  are  of  two  types; 

(1)  general  compression  of  the  building  followed  by  a  rarefaction, 

(2)  racking  of  the  building  first  in  the  direction  of  travel  of  the  boom 
and  then  in  the  reverse  direction.  The  overall  compression  of  the  build¬ 
ing  causes  an  increase  in  the  pressure  inside  tbe  building.  This  inter¬ 
nal  pressure  is  usually  one-fourtb  to  one-half  the  peek  boom  pressure  end 
significantly  sdds  to  the  structural  stiffness  in  compression. 

Window  deflections  during  a  sonic  boom  ere  e  function  of  the  exter¬ 
nal  pressure,  structure  motion,  and  internsl  pressure  built  up  through 
compression  of  the  building.  The  external  pressure  shows  evidence  of  re¬ 
flections  from  ground  and  buildings  and  pssssge  sround  buildings  and  so 
may  not  correspond  closely  to  the  free-field  boom  pressure.  To  cslculete 
the  Internsl  pressures  snd  building  mol  ^na  it  will  be  neceaaery  to  know 
(l)  motion  of  the  walls  and  roof  caused  dv  a  boom  pressure,  snd  (2)  metier, 
of  the  w*lls,  ceilljig,  and  roof  caused  by  intemel  pressure.  If  internal 
end  externel  pressure  ere  known,  window  motion  can  be  calculated  ea  the 
motion  of  a  nonlinear,  multimodal  system  under  the  action  of  known  forces. 
Such  an  analysis  would  be  a  major  effort.  Instead,  s  on e -degree -of - 
freedom  analysis  was  made  end  the  contribution  of  nonllnearltles  and 
higher  modes  computed  as  perturbations  of  the  single  mode  solution. 


7 


i 

1 

1 


IN  FREE  FIELD 


TIME 


ON  TRAILING  EDGE 


DIRECTION  OF  BOOM 


□_ 

LEADING  EDGE 

TRAILING  EDGE 


ri-nun 


FIG.  2  PRESSURES  ON  AND  WITHIN  A  BUILDING 


To  make  tbe  problem  tractable  the  following  assumptions  were  intro¬ 
duced; 

1.  The  outside  boom  Is  an  Ideal  symmetric  N-wsve,  This  assumption 
eliminates  the  rise  time  and  minor  irregularities  in  the  wsve 
form.  If  the  period  of  the  structure  is  much  longer  than  the 
rise  time  or  durations  of  Irregularities,  then  this  assumption 
Is  Justified. 

2.  The  inside  pressure  hss  the  form  of  s  full  sine  wsve  with  the 
same  duration  ss  the  boom  signature.  This  pressure  history  is 
s  simplification  of  observed  Internal  pressure  histories. 
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The  variation  In  srrlvel  times  around  the  building  cen  be 
neglected. 

The  several  ens lyses  required  ere 

•  Llneer  enelyels  of  first  mode  motion  under  sn  N-wave 

•  Llneer  analysis  of  first  mode  motion  under  Interns!  pressure 

•  Linesr  analysis  of  first  mode  motion  under  sn  N-wsve  end 

Internal  pressure 

•  Llneer  enelysls  of  multimodel  motion  under  en  N-weve 

•  Large  deflection  enelysls  of  first  mode  motion  under  en  N-wsve 

•  Calculation  of  the  amplitude  of  lnternel  pressure 

The  results  of  these  enelyses  ere  presented  next.  Following  them, 
the  results  are  collected  to  produce  sn  analysis  of  e  window  motxon  end 
stress . 

B.  Window  Deflection  and  Strees  Based  on  Llneer  Analyses 

Figure  3  Is  sn  Illustration  of  the  probeble  history  of  central  de¬ 
flection  of  s  window  (fundamental  mode  only)  to  en  N-shaped  pressure 
wave.  The  epplled  pressure  Is  shown  es  the  dotted  line  with  sn  emplltude 
corresponding  to  the  stress  which  would  occur  If  the  pressure  were  epplled 
statically.  The  dynemlc  deflection  curve  shows  oscillation  st  tbe  funda¬ 
mental  frequency  of  the  piste.  The  piste  is  oscillating  ebout  the  ststlc 
deflection  curve  so  thst  the  dynemlc  deflection  eppeers  ss  s  superposition 
of  stetlc  deflection  end  free  vlbrstlons.  In  this  figure  the  maximum 
negative  deflection  Is  approximately  equel  to  the  maximum  positive.  Tbe 
two  maxima  sre  26$  snd  40$  higher  then  the  stetlc  stress  corresponding 
to  the  maximum  spplled  pressure.  The  ratio  of  maximum  dynemlc  stress  or 
deflection  to  the  ststlc  velue  Is  referred  to  as  the  dynamic  amplification 
factor . 

The  dynemlc  amplification  fector  for  the  fundemontel  mode  of  s  equere 
plate  (applicable  to  both  stress  end  deflection)  Is  shown  In  Fig.  4.  The 
loading  wes  s  symmetric  N-sheped  pressure  weve.  In  this  figure  there  ere 
separate  curves  for  the  maximum  positive  end  negative  deflection  during 
the  boom  end  the  magnitude  of  the  free  oscllletion  following  the  boom. 
These  curves  were  derived  from  the  snslysls  In  Appendix  A  for  linesr 
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FIG.  3  PLATE  DEFLECTION  UNDER  SONIC  BOOM  LOADING 
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NONOIMENSIONAL  BOOM  OURATION.  wt/2w 

4  ri-uH-i 

FIG.  4  DYNAMIC  AMPLIFICATION  FACTOR  FOR  A  SQUARE  PLATE: 
I  LINEAR,  FIRST  MODE  DEFLECTION  UNDER  BOOMS 
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amplification  factors.  The  abscissa  is  urr/2TT,  a  nondimensions  1  relation 
between  uj,  the  fundamental  circulsr  frequency  of  the  window,  and  t,  the 
duration  of  the  N-wave. 

Higher  modes  slso  contribute  to  the  deflection  and  stress  in  the 
window.  Dynamic  smplif ication  factors  for  the  deflection  snd  streas 
under  s  symmetric  N-shaped  pressure  wsve  for  the  first  eleven  modes  sre 
shown  in  Figs,  5  snd  6.  The  first  eleven  modes  were  calculated  but  the 
msin  contributors  were  the  first  (l,l)  and  second  symmetric  (l,3)  modes. 
For  comparison  the  two  figures  include  the  positive  snd  free  vibration 
maxima  for  the  fundamental  mode.  We  note  that,  while  the  smplif ication 
factor  was  the  same  for  stress  snd  deflection  in  the  fundsmentsl  mode, 
it  differs  for  the  multimodal  case.  The  factor. for  deflection  is  modi¬ 
fied  only  very  slightly  by  the  contributions  of  higher  modes:  the  first 
pesk  is  about  4$  higher  snd  other  curves  are  modified  less.  The  ampli¬ 
fication  fsctor  for  stress  is  increased  0.60  in  some  regions  of  the  sb- 
scissa  and  the  peaks  of  the  curve  sre  broadened.  In  addition  the  second 
symmetric  mode  is  sufficiently  important  to  provide  well-defined  humps 
cm  some  of  the  curves.  On  the  sversge  the  stress  increase  during  free 
vibration  is  0.28  snd  for  the  positive  maximum  the  increase  is  0.26. 

C.  Deflection  Under  Interns!  Pressure  Based  on  Linear  Analyses 

The  response  of  a  one -degree -of -freedom  systora  to  s  pressure  in  the 
form  of  a  full  sine  wave  is  shown  in  Fig.  7.  This  is  the  expected  form 
of  internal  pressure.  Agsin  the  positive  maximum,  negative  maximum,  snd 
free  vibration  maximum  sre  shown  as  a  function  of  fundamental  frequency. 
The  amplification  fsctor  is  lsrger  than  2.0  only  near  coincidence  of  the 
forcing  frequency  snd  the  nstursl  frequency.  The  analysis  on  which  these 
curves  sre  bssed  is  detailed  in  Appendix  B. 

In  Appendix  B  sn  snslysis  is  also  msde  of  the  response  of  a  one- 
degree  -of  -freedom  system  to  a  combinstion  of  the  N-wave  snd  internal 
pressure.  Some  results  sre  shown  in  Figs.  8  and  9  for  deflection  (stress 
curves  sre  ident^csl).  The  main  effect  of  increasing  the  internal  pres¬ 
sure  is  to  greatly  decrease  the  smplif lestion  factors  in  the  vicinity  of 
«rr/2n  *  l.o. 
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NONDIMENSIONAL  BOOM  DURATION,  oir/2 v 
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FIG.  5  LINEAR,  MULTIMODAL  DYNAMIC  AMPLIFICATION  FACTOR 
FOR  CENTRAL  DEFLECTION  OF  A  SQUARE  PLATE 
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MAXIMUM  CENTWAL  MOMENT  /  STATIC  MOMENT 


NONDIMENSIONAL  BOOM  DURATION,  vr/Zir 

|4«H 

FlU  6  LINEAR,  MULTIMODAL  DYNAMIC  AMPLIFICATION  FACTOR 
FOR  STRESS  IN  THE  CENTER  OF  A  SQUARE  PLATE 
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MAXIMUM  DYNAMIC  RESPONSE /' STAT  1C  RESPONSE 


NON  DIMENSIONAL  DURATION  OF  LOADING,  ut/ 2r 


FIG.  7  DYNAMIC  AMPLIFICATION  FACTOR  FOR  FIRST  MODE  DEFLECTION 
OF  A  SQUARE  PLATE  TO  INTERNAL  PRESSURE  ONLY 
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MAXIMUM  DYNAMIC  RESPONSE  /  STATIC  RESPONSE 


MAXIMUM  OYNAMIC  RESPONSE  /  STATIC  RESPONSE 


NONDIMENFIONAL  DURATION  OF  LOADING,  wt/2tt 


'•••MM2 

FIU.  9  DYNAMIC  AMPLIFICATION  FACTOR  FOR  FIRST  MODE  DEFLECTION 
OF  A  SQUARE  PLATE  TO  BOOM  AND  INTERNAL  PRESSURE:  q2  -  0.50 
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D.  Nonlinear  Deflections  and  Stressea  Under  Booms 


In  reaponse  to  normal  wind  presaure  loadings,  window  deflection  is 
a  markedly  nonlinear  function  of  the  pressure.  This  nonlinear  behavior 
occurs  because  the  plnte  resists  part  of  the  load  as  a  membrane  after  the 
deflection  becomes  large.  According  to  the  discussion  of  Freynik  (1963) 
thia  nonlinearity  caus es  the  motion  to  differ  considerably  from  that 
shown  in  Fig.  3.  Ins  ead  of  the  smooth  sinusoidal  vibrations,  the  motion 
ahowa  abrupt  changes  u<  direction.  The  natural  frequencies  of  a  window 
undergoing  large  deflections  are  a  function  of  the  deflection  amplitude 
and  are  almost  a  continuous  band  of  frequencies  above  the  fundamental. 
With  auch  a  mixture  of  frequencies  in  the  motion,  the  natural  frequencies 
cannot  be  readily  determined  from  the  records.  The  nonlinearities  have 
an  important  effect  on 

•  atreaaea  and  deflections, 

•  natural  frequency, 

•  dynamic  amplification  factor. 

We  will  consider  each  of  theae  effects  and  then  provide  some  numerical 
values  to  guide  in  aaaesalng  the  seriousneas  of  the  nonlinearity. 

The  variation  of  streaa  and  deflection  with  applied  static  pressure 

ia  ahown  in  Fig.  10  baaed  cxi  the  calculations  in  Appendices  C  and  D.  In 

Appendix  C  the  forms  of  the  equations  are  derived  theoretically.  The 
coefficients  of  the  equations  are  evaluated  in  Appendix  D  from  experimen¬ 
tal  data  of  other  ..investigations.  Figure  10  ahows  the  nondlmenslonal 
quantities  w0/h  (deflection,  ^),  oa2/Eha  (stress,  S),  and  qa4/Eh4  (pres¬ 
aure,  Q),  where 

w0  =  central  deflection  of  the  square  plate 
a  =  stress  (bending,  membrane,  or  both) 

q  =  uniform  pressure  on  the  plate 

a  =  side  of  plate 

h  =  thickneaa  of  plate 

E  =  elastic  modulus  of  plate  material 
The  nonllnearltlea  became  important  for  a  deflection  greater  than  0.5 
times  the  plate  thickneaa.  A  2  paf  boom  will  cauae  deflections  up  to  1.9 
times  the  thickneaa  in  aorae  windows  -  that  ia,  well  into  the  nonlinear 
range . 
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NONOIMENSIONAL  STRESS, 


NONOIMENSIONAL  PRESSURE,  Q  ■  <jo4/Eh4 


FIG.  10  DEFLECTION  AND  STRESS  IN  THE  CENTER  OF  A  SQUARE  PLATE 
AS  A  FUNCTION  OF  STATIC  PRESSURE 


In  Fig.  11  it  ia  apparent  that  tha  relation  between  deflection  and 
central  stress  is  not  altered  much  by  the  nonlinearities .  The  fact  that 
the  linear  and  nonlinear  curves  do  not  coincide  ct  the  origin  suggests 
that  there  may  be  some  uncertainties  about  the  experimental  values  from 
Appendix  D  on  which  the  curves  are  based. 


nonoimensicnal  deflection,  {  *  w0/h 

TS-60M-7 

FIG.  II  NONLINEAR  RELATION  BETWEEN  CENTRAL  DEFLECTION  AND  STRESS 


Next  we  consider  the  large  deflection  of  a  window  under  boom  pres- 
aures.  Dynamic  calculations  of  plates  undergoing  large  deflections 
normally  begin  with  von  Karman'a  equations.  (See  Timoahenko  (1936)  for 
a  derivation  of  the  equations.)  However,  the  reaulta  from  von  Karman'a 
equations  did  not  agree  well  with  experimental  data  on  window  deflections 
(aee  Appendix  D).  Therefore  we  took  the  preaaure-def lection  relation 
derived  froo  the  static  experiments,  restricting  attention  to  motion  in 
the  fundamental  mode  only. 
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Nonlinear  deflection  of  a  window  under  an  N-wave  Is  shown  In  Figs. 

12,  13,  and  14  (based  on  calculations  In  Appendix  E).  These  three  figures 
correspond  to  peek  pressure  amplitudes  required  to  cause  static  deflec¬ 
tions  of  h,  2h,  end  3h.  (£  Is  the  nondlmenslonal  static  deflection 

ws/h  under  the  peek  pressure.)  The  peaks  of  these  figures  and  Fig,  4  are 
summarized  In  Fig.  15.  In  the  summary  figure  we  see  that  large  deflec¬ 
tions  tend  to  reduce  the  period  of  the  plate  ao  that  the  maxima  occur  at 
smaller  values  of  u)t/2tt  than  occur  from  small  deflections.  Also,  the 
amplification  factor  for  free  vibration  Is  reduced  for  large  deflections 
but  Is  unaffected  for  the  positive  maximum  during  forcing.  The  change 
In  the  apparent  natural  frequency  with  deflection  Is  summarized  In  Fig.  16. 

An  estimate  of  the  multimodal  nonlinear  behavior  Is  made  In  Appendix 
E  for  both  stress  and  deflection.  The  peak  deflections  are  modified  very 
slightly  by  the  pertlclpatlon  of  modes  above  the  fundamental.  Therefore, 

It  Is  recommended  that  the  nonllneer  response  In  the  fundamental  mode  be 
used  to  calculate  deflections  with  no  changes  to  account  for  higher  modes. 


The  estimate  for  peak  stress  Is  more  complicated,  and  the  accuracy 
of  the  estimate  may  be  considerably  Jess  than  that  for  deflection.  The 
estimation  procedure  to  be  outlined  eppeers  loglcel  but  has  not  been 
verified  In  any  way.  The  procedure  Includes  two  elements; 

•  Conversion  of  the  nonlinear  deflection  amplification  factor 
to  stress  amplification  factor 

*  Addition  ol  an  increment  to  correspond  with  the  contribution 
of  higher  modes. 


To  convert  from  deflection  to  stress,  Eq.  D.8  can  be  used  to  relate 
the  stress  In  the  first  mode,  Sx,  to  the  totel  central  deflection, 

S  .  =  4  pr  M  4.  fl  1K7l!l 

A  \  -  -• - *3/ 

But  the  deflection  E  Is  F_?  where 

?  s 

Fp  Is  the  amplification  factor  for  deflection  and 
E  Is  the  static  deflection. 

Therefore 

Sj  =  4.9  F  5  (l  +  0  167  FA  ) 

'  8  5  s 

The  next  step  Is  to  add  a  f  utor  corresponding  to  the  contribution  of 
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NONOIMENSIONAL  BOOM  DURATION.  wr/Zw 


FIG.  12  NONLINEAR  FIRST  MODE  DEFLECTION  OF  A  SQUARE  PLATE 
TO  A  BOOM  LOADING  WITH  AMPLITUDE  SUCH  THAT  £s  -  1.0 
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NONOIMKNSIONAL  BOOM  OURATION,  iht/Zv 


f!G.  13  NONLINEAR  FIRST  MODE  DEFLECTION  OF  A  SQUARE  PLATE 
TO  A  BOOM  LOADING  WITH  AMPLITUDE  SUCH  THAT  -  2.0 
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NONDIMENSIONiL  BDOM  DURATION,  ur/2 w 
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FIG.  15  SUMMARY  OF  NONLINEAR  FIRST  MODE  DEFLECTION 
OF  A  SQUARE  PLATE  TO  A  BOOM  LOADING 


0  0.5  1.0 

_ NATURAL  PERIOD _ 

NATURAL  PERIOD  FOR  ZERO  DEFLECTION 

TA-*0*S-I2 


FIG.  16  VARIATION  OF  APPARENT  NATURAL  FREQUENCY 
WITH  CENTRAL  DEFLECTION 


the  higher  modes.  According  to  the  sns iysis  of  Appendix  E,  the  stress 
sssocisted  with  these  modes  is  primsrily  bending  (not  membrsne)  stress 
and  hence  equsl  to  thst  obtsined  for  the  linesr  case.  Therefore  sdd  0.28 
times  the  totsl  linesr  static  stress  if  the  peak  occurs  during  free  vi- 
brstion,  or  0.26  if  the  pesk  occurs  during  forcing.  That  is,  for  free 
vibration  the  dynamic  stress  is 

S  *  4.9  F.|  (1  +  C.167  F.S;  )  +  0.28  (4.4)  $ 

S  s  5,  s  3 

The  stress  smpllf icstion  factor,  F  ,  is  found  by  dividing  S  by  S  where 

cj  £ 

S  =  4.4  5  (1  +  0.186  F  ) 

a  Ss  'S 
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Aa  an  example  conalder  a  aquare  window  and  a  boom  with  the  follow¬ 
ing  character iatlca 

Window  area  =  45  aquare  feet 

Thlckneaa  =  l/4  Inch 

Natural  frequency  =  7 .4  Hz 

Boom  duration  =»  0.1  seconds 

Peak  pressure  -  8.85  paf 

Nondimensional  pressure  =  66 

According  to  Fig.  10,  the  atatlc  deflection  under  the  peak  preasure  la 
2h.  F^j  is  1,80  according  to  Fig.  15,  using  an  abaclaaa  of  0.1  (7.4)  = 

0.74,  Then  for  moment  we  compute 

4.9  (1.80)  (1  +  0.167  x  2.0  x  1.80)  +  0.28  (4.4) 
a  '  4.4  (1  +"0.186  x  2.0) 

=  2.54  (nondlmenalonal) 


a  value  which  appeara  quite  reaaonable. 


Now  let  us  take  three  examples  of  wlndowa  and  determine  the  deflec¬ 
tions,  stresaea,  and  frequencies  In  response  to  a  2  psf  boom,  disregard¬ 
ing  internal  preasure.  The  maximum  and  minimum  slenderneaa  ratios  (340 


and  220)  and  an  intermediate  value  will  be  used  for  the  comparison.  The 
calculated  values  are  given  in  Table  3.  Included  in  the  table  are  values 
from  the  linear  analysis  of  the  same  problem.  The  window  with  the  largest 
a/h  value  gives  the  most  markedly  nonlinear  behavior.  For  thia  window 
the  large  deflection  theory  gives  a  deflection  juat  84$  of  that  from  the 
linear  theory,  and  stresses  75$  of  those  from  the  linear  theory.  The 


windows  with  smaller  a/h  ratios  give  deflections  about  equal  to  thr* 


linear  values  although  the  stresses  are  considerably  lower. 
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TABLE  3 


DEFLECTION  AND  STRESS  OF  SQUARE  WINDOWS  UNDER  A  2  PSF  BOOM 


a/h 

Static  Loading 

Dynamic  Loading 

Linearized  Analysis, 
Dynamic  Loading 

w0/h 

Deflection 

a 

Stress 

(psi) 

w0/h 

Deflection 

a 

Stress 

(psi) 

w0/h 

Deflection 

a 

Stress 

(pai) 

220 

0.15 

160 

0.30 

370 

0.30 

440 

280 

0.39 

240 

0.75 

560 

0.78 

730 

340 

0.78 

350 

1.43 

810 

1.70 

1080 

NOTE:  The  boom  duration  and  fundamental  frequency  of  the  window 

were  preaumed  to  be  auch  that  ;ht/2tt  =  1.0. 


E .  Determination  of  the  Amplitude  of  Internal  Preasure 

The  preaaure  In  a  building  which  la  struck  by  a  boom  may  be  caused 
by 

1.  flow  of  the  preasure  weve  through  openings  In  the  building 
{open  doors,  windows,  ventilation  ports), 

2.  overall  compression  of  the  building, 

3.  transmission  through  flexible  areas  such  as  windows. 

The  third  Is  essentially  a  different  way  of  stating  the  second  with  the 
added  presumption  that  the  only  elements  which  will  deflect  appreciably 
are  windows.  The  third  cause  was  assumed  by  Blume  (1965)  to  be  the  most 
important.  In  this  section  a  calculation  is  presented  determining  the 
pressure  rise  as  a  function  of  overall  building  compression.  The 
calculation  showb  that  building  compression  may  be  the  prims  cause  of 
Internal  pressure. 

For  the  calculation  we  will  assume  that  the  history  of  the  internal 
pressure  la  sinusoidal  with  period  equal  to  the  boom  duration.  Internal 
pressure  is  calculated  from  the  equilibrium  condition  between  the  boom 
preaaure,  wall  deflection  ond  internal  pressure.  The  necessary  equations 
are  derived  using  the  variables  illustrated  in  Fig.  17.  The  maximum 
deflection  of  the  roof  is 
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(1) 


FIG.  17  PRESSURES  CONSIDERED  IN  CALCULATING 
THE  INTERNAL  PRESSURE  OF  A  BUILDING 
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Comparable  equations  can  be  written  for  the  wall  and  window  deflections. 
Similarly,  the  celling  deflection  Is 


wc  =  Fic  V  (q’>  '  q’2) 


(2) 


where 

Fjc  Is  tl»e  dynamic  amplification  factor  for  the  celling  under 
internal  pressure 

w0„  is  the  average  static  celling  deflection  under  unit  pressure 
q23  la  the  pressure  In  the  first  floor  rooms. 

The  internal  pressure  Is  caused  by  deflection  of  the  structural  surfaces 
so  that 

-  wc°c  *  Vw  *  WB°B  (3) 

q»*  =  Vr  -  WCGRC  (4) 

where 

wtf  =  window  deflection 
w^  wall  deflection 

G  ,  Gtf,  GQ,  Gp,  and  GRC  are  factors  relating  the  average  deflection 
to  the  change  In  pressure. 

The  main  contribution  to  internal  pressure  Is  from  large  panel  elements 
such  as  the  roof,  large  doors,  or  windows. 

To  validate  the  method,  a  calculation  waa  made  of  the  peak  Internal 
pressure  Jn  a  simple  rectangular  one-atory  structure  designated  PF-6  In 
the  White  Sands  sonic  boom  teat  report  of  John  A.  Blume  and  Associates 
(1965). 

The  experimental  results  obtained  during  the  sonic  boom  tests  Indi¬ 
cate  a  nonlinear  relation  between  Internal  and  external  pressure.  For 
low  presaures  the  ratio  Is  about  40$  for  booms  from  an  F-104,  At  5  psf 
the  ratio  Is  30$.  For  our  calculation  the  2  x  4's  in  the  walls  and  roof 
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were  aoaumed  to  provide  all  the  atlffnera  to  those  elements  and  to  be 
fixed  at  both  ends.  Yhe  poaitlve  maximum  amplification  factor  waa  used 
from  Figs.  4,  7-9.  The  boom  duration  was  assumed  to  be  0.10  seconds,  a 
nomint  1  value  for  the  F-104 .  Calculated  frequencies  of  the  atructural 
elements  were 

Windows  18  Hz 

Walla  25  Hz 

Roof,  Celling  8  Hz 

Equations  1  through  4  were  solved  for  a  1  paf  boom  to  give 
q21  =  0.64  psf  (peak  attic  pressure) 

q22  =  0.36  psf  (peak  pressure  In  first  floor  rooms) 

The  value  of  0.36  for  q22  should  be  compared  to  the  Internal  preaaure  of 
30  to  40$  measured  by  Blume  and  Associates, 

The  correspondence  between  the  calculated  and  experimental  valuea  la 
very  good  considering  the  uncertain  basis  for  aome  of  the  assumptions. 

The  present  calculations  alro  ahow  the  relative  importance  of  the  contri¬ 
bution  of  each  element  to  the  pressure  rise  In  the  structure.  For  PF-6, 
the  pressure  caused  by  motion  of  the  roof  and  celling  waa  82$  of  the  In¬ 
ternal  pressure,  q22.  The  walls  contributed  8$,  windows  10$.  While  the 
calculations  are  not  verified  sufficiently  to  use  for  predictions,  the 
results  led  us  to  believe  that  compression  Is  the  main  cause  of  the  rise 
In  internal  pressure. 

The  preceding  procedure  is  dependent  on  the  validity  of  certain  as¬ 
sumptions  concerning  *>nd  fixity  of  wall  and  roof  members,  and  on  win¬ 
dows  and  doors  being  closed.  Instead  of  conducting  this  procedure,  which 
Is  unverified,  it  Is  reasonable  to  depend  on  an  estimate  of  the  internal 
pressure  based  on  the  experimental  results  of  Blume  (1965).  For  a  flexi¬ 
ble  structure  the  internal  pressure  might  be  estimated  at  50$  of  the  boom 
pressure.  If  the  roof  or  celling  or  upper  floor  la  a  stiff  element  and 
the  walls  are  not  abnormally  flexible,  an  educated  guesa  lndlcatea  that 
the  Internal  presaure  is  25$  or  leas  of  the  boom  pressure.  Following 
the  reduction  of  the  Edwards  Air  Force  Base  test  data  It  la  expected  that 
a  more  reliable  baais  for  these  pressure  estimates  will  be  available. 
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F .  Summary  of  Window  Behavior  Under  Boom  Loading 


Determination  of  the  response  of  windows  to  booms  requires  the  com¬ 
putation  of  a  host  of  factors.  Now  that  each  of  the  factors  has  been 
introduced,  a  summary  of  the  calculation  is  presented: 

1.  Eatimate  window  deflection  using  some  guess  at  internal  pres¬ 
sure;  use  Figs.  4,  8,  and  9.  Thia  estimate  is  required  for 

Step  2. 

2.  Compute  the  internal  pressure  acting  on  the  inside  of  the  win¬ 
dow  as  outlined  in  Section  E. 

3.  Determine  the  amplification  factor  for  the  window  as  a  one- 
degree  -of  -freedon  syatem  to  boom  pressure  and  internal  pressure 
use  Figs.  4,  8,  and  9. 

4.  Modify  the  computed  factor  for  nonlinear  effects:  Find  the 

dynamic  amplification  factors  for  appropriate  static  deflection 
and  for  zero  deflection;  multiply  the  amplification  factor  from 
Step  3  by  the  ratio  of  these  factora  to  obtain  Fr  . 

5.  Modify  the  values  for  the  multimodal  effect:  For  deflection, 

F. 


let  F^  = 


5»' 


For  moment  or  stress,  compute  F^;  the  dynamic 


amplification  factor  fi*om  the  relation 


F 


1.115  Fff  (1  +  0.167  Fet  )  +  0.28 

_ 5 _ L® _ 


1  +  0.186  ? 

3 

The  valuea  of  streaa  and  deflection  can  then  be  found  from  the  equations 


5  =  F  A 

5  s 


and 


S  =  FS  =  4.4  E  ( 1  +  0 . 186  E  )  F 

os  5s  Ss  a 

It  ia  expected  that  the  dominant  displacement  response  of  square 
windows  to  a  2  psf  far -field  sonic  boom  will  be  resonance  in  the  funda¬ 
mental  frequency  of  the  window.  For  longer  rectangular  windows  there 
will  be  oacillation  in  the  first  and  second  symmetric  modes.  The  noar- 
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field  boom  exhibits  oscillations  which  may  excite  oth~r  window  frequencies. 

Nonlinear  effects  are  expected  for  windows  with  large  slenderness 
ratios  (length  to  thickness  over  300)  when  subjected  to  2  psf  booms. 

Instead  of  a  resonance  in  the  fundamental  mode  at  the  fundamental  fre¬ 
quency,  the  resonance  will  occur  at  a  range  of  frequencies  near  the  fun¬ 
damental.  Freynik  (1963)  has  shown  that  the  Fourier  spectrum  of  response 
has  a  continuous  band  of  frequencies  with  the  frequencies  from  linear 
theory  dominant .  The  linear  frequencies  hecome  less  dominant  as  the 
deflections  increase.  The  magnitude  of  the  peak  deflection  and  stress 
will  be  reduced  by  the  nonlinearities.  The  dynamic  amplification  factor 
for  deflection  and  stress  will  be  1.4  to  2,7. 

Maximum  stress  and  maximum  deflection  are  not  closely  related,  even 
in  the  static  case.  Maximum  deflection  under  a  uniform  static  pressure 
is  essentially  a  function  of  the  first  mode  response  whereas  maximum 
stress  receives  a  contribution  of  about  10^,  from  the  higher  modes.  In 
the  case  of  dynamic  and  particularly  nonlinear  dynamic  response,  there 
is  even  less  correlation  between  the  stress  and  deflection.  Hence,  if 
stress  is  desired  from  experimental  data,  stress  should  be  determined 
from  strain  gage  measurements,  and  if  deflection  is  desired,  deflection 
t.iould  be  measured.  One  should  not  depend  on  calculating  one  quantity 
from  a  measurement  of  the  other  at  a  given  point. 
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SECTION  IV 


STATISTICAL  PREDICTION  OF  FAILURE 


Tha  most  Important  information  which  should  be  obtained  from  a  study 
of  window  raaponsa  to  aonic  booms  ia  an  answer  to  the  question:  What  is 
tha  probability  of  failure  of  win  ,  ’*»  undar  aonic  boom  presaurea?  To 
anawar  this  quaation  wa  naad  ataxii;;  leal  data  of  three  types: 

•  Variability  of  raaasured  .-.xlmum  prasaures  for  a  planned  pres¬ 
sure  of  the  sonic  boom.  Thesa  variations  are  dependent  on 
weather  conditions,  position,  velocity,  and  acceleration  of 
aircraft,  and  terrain.  With  our  present  state  of  knowladga 
about  booms,  these  variations  muat  be  handlad  statistically 
rather  than  deterministically. 

•  Distribution  of  window  sizas  in  buildings. 

•  Distribution  of  window  strengths  in  normal  kinds  of  caaings. 

We  hava  made  two  computations,  aaauming  flrat  that  window  glasa 
strangtha  are  distributed  normally,  and  second  that  they  are  distributed 
lognormally,  (Available  data  do  not  distinguish  between  theae  distribu¬ 
tions.)  We  take  the  coefficient  of  variation  of  the  strength  as  25$  as 
recommended  by  McKinley  (1964)  of  Pittsburgh  Plate  Glasa.  The  mean  boom 
strength  is  taken  as  2  psf  with  a  coefficient  of  variation  of  25$.  The 
25$  is  characteristic  of  fairly  calm  daya .  The  paak  stresses  used  will 
ba  thoaa  from  Tabla  3,  for  a/h  =  340  (moat  critical  case)  and  for  a  glaaa 
strength  of  6600  pai . 


The  probability  of  failure  is  treated  from  the  viewpoint  of  multiple¬ 
valued  random  phenomena  as  described  by  Parzen  (1960)  (Chap.  7,  Sec.  3): 


<p2(x2j  dx2  dXj 


where  cp's  are  mutually  independent  probability  density  functions.  For  a 
normal  distribution,  9  is  given  by 
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and  for  the  lognormal,  cp  is  (see  Aitchison  and  Brown  (1957)) 


1 

cp(x)  =  / 2ttxc  exp  [-1/2  (l0gQX  m)  2j 

While  an  exact  analytical  solution  of  Eq.  5  ia  not  possible,  the  integra¬ 
tion  can  be  easily  performed  numerically.  These  computations  were  made 
for  the  moat  critical  caae  of  Table  3,  that  ia,  for  a/h  -  340.  The 
results  are 


Distribution  Probability 

Normal  .0002 

Lognormal  10-9 

Evidently  then  the  failure  probability  is  critically  dependent  on  the 
distribution  aasumed.  It  is  often  aaid  that  the  normal  distribution  is 
impoaaible  because  it  predicta  a  finite  poasibility  of  negative  values 
of  atrength  and  boom  pressure.  However,  if  we  truncate  the  normal  dis¬ 
tribution  at  zero  for  the  above  calculation,  we  make  no  noticeable  change 
in  the  calculated  probability.  Hence  there  appears  to  be  neither  logical 
nor  experimental  baaea  for  determining  the  correct  distribution. 

From  claima  data  it  appears  that  the  probability  of  damage  for  win- 
dowd  per  2  paf  boom  is  of  the  order  of  10_a.  This  indicates  that  the 
correct  diatribution  for  window  strengths  may  be  intermediate  between  the 
normal  and  lognormal.  Because  of  the  very  low  probability  value  of  10“6, 
it  ia  not  to  be  expected  that  laboratory  tests  can  provide  the  distribu¬ 
tion  with  aufficient  accuracy  for  damage  calculations. 
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SECTION  V 
SUMMARY 


The  results  of  this  report  apply  only  to  windows  which  may  be  con¬ 
sidered  simply  supported.  An  unknown  portion  of  reel  windows  sre  looee 
or  have  stress  reiaers  in  their  mountings.  Most  windows  without  such 
defects  will  feil  under  boom  pressures  of  100  psf  or  less  while  st  least 
e  few  will  fsil  et  10  or  20  psf.  Defective  windows  mey  fail  at  even  low¬ 
er  levels.  At  stetic  pressures  of  1  to  10  psf  the  response  of  simply 
supported  windows  becomes  nonlineer.  The  nonlineerity  radically  effects 
the  reletions  emong  deflection,  strees,  and  pressure. 

The  dominent  motion  of  s  window  under  e  2  psf  boom  ie  oaoilletion  in 
the  fundamental  frequency.  For  windowe  with  s  length  to  thickness  retie 
over  300,  nonlineer  effects  are  expected  to  increese  the  response  in  the 
higher  modes.  Also  the  motion  in  the  fundemental  mode  will  occur  at  e 
renge  of  frequencies  near  the  fundementel  frequency. 

A  procedure  was  developed  for  predicting  the  deflection  and  streas 
of  windows  under  low-pressure  booms.  The  present  enalysie  ia  adequate 
elthough  approximate.  The  procedure  is  based  on  e  one-degree-of-freedom 
enelysis  plus  estimates  of  the  multimodel  snd  nonlineer  effecta  end  of 
the  interection  with  the  building  motion. 

The  celculetion  procedure  provides  predictiona  of 

•  Peek  interne  1  pressure  esused  by  incidence  of  a  boom  on  a 
building 

•  Response  of  e  one-degree-of-freedom  system  to  e  combined  boom 
snd  internsl  pressure  loading 

•  Contribution  of  higher  modes  to  the  streas  end  deflection  of 
windows 

•  Modification  of  window  response  ceused  by  large-def laction 
(nonlinear)  effecta. 
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The  possibility  of  window  failure  caused  by  2  psf  booms  is  consi¬ 
dered  from  a  statistical  standpoint.  But  no  reasonable  estimate  of  the 
statistics  of  failure  can  be  made  because  statistical  distribution  of  the 
strength  of  glass  is  not  known  precisely. 

In  experimental  studies  of  window  motion  it  is  necessary  to  meaaure 
(a)  central  deflection  of  the  window,  (b)  strains  at  aeveral  points  on 
both  sides  of  the  pane,  (c)  pressure  on  hoth  sides  of  the  window, 

(d)  edge  fixity  of  the  window.  In  dynamic  studies  the  history  of  each 
of  the  measurements  a,  b,  and  c  must  be  taken.  Auxiliary  measurements 
must  be  made  of  Young's  modulus  and  Poisson's  ratio  and  the  exact  dimen¬ 
sions  of  the  glaas  pane. 

A  basis  for  reduction  of  experimental  data  on  window  response  is 
provided.  Stress,  deflection,  and  pressure  should  be  nondimensionalized 
and  graphed  againat  each  other  to  show  the  trends.  Dynamic  results  will 
require  a  dynamic  analysia  similar  to  those  herein  to  provide  an  adequate 
evaluation  of  the  data. 

There  ia  need  for  empirical  relations  among  stress,  deflection,  and 
preaaure  for  windowe  in  the  uaual  mountings.  Such  relations  should  pro¬ 
vide  a  baais  for  multimodal  response  calculations.  The  next  step  is  a 
dynamic  multimodal  calculation.  With  static  relations  and  dynamic  cal¬ 
culations.  an  adequate  base  is  laid  for  reducing  data  from  sonic  boom 
te3ts  on  windows. 
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SECTION  VI 


FURTHER  STUDY 

To  more  clearly  define  the  response  of  windows  to  sonic  booms,  four 
further  studies  may  be  considered.  The  first  three  should  be  undertaken 
in  the  order  lasted. 

1.  An  experiments?,  determination  of  the  static  and  dynamic  load- 
deflection  relations  for  windows.  Measurements  must  be  suffi¬ 
ciently  complete  to  provide  information  on  the  behavior  in 
several  modes. 

2.  A  calculation  of  the  dynamic  multimodal  response  of  windows 
undergoing  large  deflections. 

3.  A  detailed  analysis  of  the  response  of  windows  in  test  houses 
at  White  Sands  or  Edwards  Air  Force  Base.  The  required  instru¬ 
mentation  consists  of  an  accelerometer,  strain  gages  or  dis¬ 
placement  gages  on  the  window,  and  pressure  gages  on  the  inside 
and  outside  of  the  window.  In  addition,  a  motion  gage  on  the 
window  frame  would  be  useful  to  show  the  amount  of  support 
movement . 

4.  A  study  of  the  magnitude  of  internal  pressures  in  huildlnga 
subjected  co  booms.  The  internal  pressure  should  be  measured 
and  calculated  from  the  motion  of  structural  elements.  A  by¬ 
product  of  such  a  study  would  be  an  accurate  method  to  calcu¬ 
late  the  frequencies  and  deflections  of  building  elements  in 
residential  construction. 

While  each  of  the  above  studies  has  scientific  and  engineering  merit,  none 
appears  to  be  justified  by  the  needs  of  the  sonic  boom  program. 


APPENDIX  A 


DYNAMIC  AMPLIFICATION  FACTOR:  LINEAR 

The  dynamic  amplification  factors  for  the  stress  and  deflection  at 
the  center  of  a  simply  supported  square  plate  were  calculated  for  a 
pressure  pulse  in  the  form  of  an  N-wave r  The  static  moment  and  deflec¬ 
tion  are  given  by  Timoshenko16  (1959)  in  the  following  series  form; 

00  00 

V  V  m2-Hm2  . 

L,  **  ran  ( m  !+nJ 'l  J  8in  (mTIX/a)  sin  (nTty/a)  (A.l) 

m—  1  m=l  1 

00  00 

L  £  mniV+n4)4  sin  (m71x/a)  sin  (nny/a)  (A. 2) 

^  l  in™"  X 

where  Mx  is  the  moment  in  the  x-direction  in  the  plate 
w  is  the  deflection 
q0  is  the  pressure 

a  is  the  length  of  one  side  of  the  plate 
D  is  the  flexural  rigidity  of  the  plate 
x,y  are  coordinates  of  the  plate 
v  is  Poisson's  ratio 

m,n  are  indices  denoting  the  number  of  modes  in  the  x  and  y  directions 
respectively . 

For  purposes  of  the  computer  program  required  to  evaluate  the  equa¬ 
tions,  nondimens ionalized  static  deflection  and  moment,  ZO  and  MO,  were 

w 

defined  by 


and 


(l+  v ) 


where  w0  and  M0  ere  the  central  deflection  and  moment  under  a  uniform 
static  pressure.  Therefore 
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S;(m,n)  +  S1(n,m)  = 


mnirrr+n' 


I 


ii  Ms 


Therefore,  let 


S(ra,m) 


and 


Then 


o  t 


(a,n) 


m-n 


HOZL. 

mn(m2+n2) 


for  m  =  n 


HO  =  EE  S(m,n) 
m=l  n=m 
odd  only 


7he  dynamic  deflection  is  obtained  uaing  the  mode  auperposition 
approach  which  can  be  found  in  many  texta  (aee,  for  example  Norri3,  et  al 
(1959)).  In  thia  method,  the  deflection  is  represented  as  a  seriea  of 
products  of  three  elements; 


»(t) 


Y  Y  F  (t)  cp  (x,y) 

mn  Tmn  Ymn  * 

m  n 


(A.3) 


where 

F  (t)  ia  the  dynamic  amplification  factor 
mn 

dr  is  a  participation  factor  for  each  mode 

mn 


m  ia  the  mode  shape. 

Ymn 


For  the  square  plate 


<Pmn(x’y) 


'mn 


=  sin  (mrrx/a)  sin  (niry/a) 


a  a 


<?o« 


r 


HO  HO 


9„(«.y)  dxdy 

rmn 


a  a 
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^  (ma+na) 2 


(A. 6) 


where 


g  is  the  acceleration  of  gravity 

y  is  the  weight  per  unit  area 

oi  is  the  natural  circular  frequency  of  the  mn  mode 
mn 

q0  is  the  maximum  of  the  applied  pressure 


In  the  subsequent  analyses,  we  will  develop  two  solutions:  one 
applicable  during  forcing  and  the  other  applicable  in  the  free  vibration 
following  forcing.  Let  the  loading  be  given  by  the  expression  q0f(t). 
Then  the  dynamic  amplification  factor  i3 


F  (t)  =  oi  f(t')  ain  [oi  (t-t')]  dt 1 

mn  mn  mn 

,Jo 


(A.7) 


For  an  N-wave  the  loading  is 


q0(l-2t/r)  for  0  £  t  £  T 


where  t  ia  the  duration  of  the  N-wave.  Then  the  solution  during  the  time 
of  forcing  is  given  by 

F  (t)  *  l-cos(o)  t)  -  ^  +  — -  sin  U  t)  (0  £  t  £  t)  (A. 8) 

mn  mn  t  oi  T  mn 


j,  _  16<fo°4 

’mn  ~  TTdDmn7ma+n2y2 


for  m  and  n  odd 
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w(t)  =  1^g— -  Yj  Z  ain  (mTTX/a)  siJl  (nny/a)  F^t)  (A. 10) 
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for  (0  S  t  S  t) 
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These  values  of  w(t)  and  are  valid  only  during  the  loading,  that  Is, 
up  to  t  =  T. 

Using  the  quentities  R  and  S,  and  introducing  nondimensionallzed 
central  deflection  Z  and  moment  M,  Eqs .  A. 10  and  A. 11  take  the  following 
form. 

CO  00 

Z  =  V  Y  R(m,n)  [l-cos(u)  t )  -  ~  +  — —  sin(u)  t)]  (A. 12) 

L-‘.  *  mn  '  t  u!  T  v  mn  ' 
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00  00 

M  =  Y  Y  s(mJn)  [--cos(u)  t)  -  ^  +  — - —  sin(u)  t)]  (A. 13) 
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odd  only 

Figure  3  shows  the  history  of  the  stress  or  deflection  in  the  funda¬ 
mental  mode  during  forcing  and  during  free  vibration.  Both  for  deflec¬ 
tion  and  moment  in  the  plate,  the  major  contribution  is  given  by  the 
first  mode  (m  =  n  =  l).  Using  this  fact  we  can  approximately  locate  the 
maxima  of  the  deflection  and  moment  by  finding  the  maxima  of  the  firat 
mode.  The  maxima  occur  at  a  time  when  the  derivative  of  the  temporal 
term  is  zero,  i.e.  when 


y  sin  lot  =  1  -  COS  (tut)  (a.  14) 

This  relation  is  satisfied  in  two  ways.  The  positive  maxima  (defining 
inward  deflection)  are  given  by 


u,'nt  =  2  arc  tan  (uujit/2)  +  2  lrr  (A.15) 

max 

where  i  is  zero  or  a  positive  integer.  The  case  i  =  0  ia  of  interest 
here  because  it  defines  the  largest  of  the  positive  maxima.  The  negative 
maxima  (defining  outward  deflection)  occur  for  l^max  =  w^ere 
is  a  positive  integer.  To  obtain  the  largest  negative  value,  i  is 
chosen  so  that 

2  in  <  wT  <  2(1+1)  n  (A. 16) 
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The  maximum  values  obtained  from  Eqa.  A. 12  and  A. 13  represent  the 
pv_k  of  the  first  inward  motion  of  the  window.  To  reduce  these  values  to 
the  dynamic  amplification  factors  we  divide  them  by  the  corresponding 
static  valuea  from  Eqa.  Al.  and  A. 2.  That  ia 


A  or  C 


and 


F 

ra 


MO 


B  or  E 


where  F^  and  F^  stand  for  dynamic  amplification  factors  for  deflection 
and  moment.  A,  B,  C,  E,  are  names  for  theae  factors  in  the  computer 
program  developed  to  evaluate  the  equations.  A  and  B  are  for  positive 
maxima,  c  and  E  for  negative  maxima. 

In  addition  to  finding  the  amplification  factor  for  the  multimodal 
case,  it  ia  of  intereat  to  find  the  factor  for  the  firat  mode.  This 
calculation  is  eaally  made  becauae  the  Eqa.  A.l^  and  A. 16  exactly  locate 
the  maxima.  The  amplification  factora  for  deflection  and  moment  in  the 
firat  mode  only  are  dealgnated  Al  and  Cl.  Curvea  showing  the  amplifica¬ 
tion  factora  for  deflection  and  moment  are  ahown  in  Figs.  4,  5,  and  6. 

Following  the  application  of  the  load  the  plate  undergoes  free  vi¬ 
brations  which  are  also  associated  with  deflections  larger  than  the 
static  onea .  The  deflections  during  a  period  of  free  vibrations  have 
the  form 


w  =  w0  coa  U(t-t0)]  +  sin  [uu(t-t0)] 


<JU 


...  ~  4  —  *  *■  4  -  T  J  «  4>  1  4  J  . J*  1 .4  4..  -4.  it-  4. 

»»  Q  diiU  *1  Q  Oi  O  AUJLt  lOl  UCllOLU^Ii  ClfrlVJ  VCiUbl  bj  (4  1/  UiC  I'llilV  «-  Q  # 

Therefore,  the  firat  step  in  calculating  these  free  motions  is  to  deter¬ 
mine  the  displacement  and  velocity  at  the  termination  of  forcing  and 
beginning  of  free  vibration,  that  is,  at  t  =  T.  The  displacement  and 
velocity  at  the  center  of  the  plate  are 
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The  equations  for  moment  efter  the  loeding  cen  be  determined  by  the 
same  procedure.  The  maxima  in  the  first  mode  ere  given  by 

ten  Ouft-T)]  =  — y (A. 22) 

"  1  ill)!  i 

which  reduces  to 

Wilt  -  ^  +  i(n)  (A. 23) 

where  i  is  zero  or  e  poeitive  integer,  i  must  be  lerge  enough  eo  thet 

t  >  T. 
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The  dynamic  amplification  factor  for  the  firat  mode  is  deaignated 
A2  and  is  obtained  using  the  tima  given  by  Eq.  A. 23.  The  multimodal 
amplification  factors  are  designatod  APR  and  MOMFR,  where 

APR  -  ^  end  MOMFR  =  ^ 

Z  and  M  are  the  maximum  values  of  nondimensional  deflaction  and  moment 
which  occur  during  free  vibration. 

Flow  charts  for  the  program  are  shown  in  Figs.  A.l  -  A. 3.  Most  of 
the  variables  have  been  defined  in  the  previous  discussion.  Comment 
carda  in  the  program  listing  further  help  to  clarify  the  computation 
ayatem  uaed.  The  program  ia  written  in  Fortran  IV  and  has  been  uaed 
with  both  a  Burrougha  B5500  and  a  CDC3200.  The  input  parameters  required 
for  the  program  are  explained  in  the  program.  A  sample  set  of  data  cards 
are 

(1)  2 

(2)  40  21  0.02  0.02  0.005 

(3)  20  21  1.00  0.05  0.005 

The  firat  card  indicates  that  there  will  be  two  groupa  of  data.  The 
second  requires  a  calculation  for  40  values  of  the  abaciaaa,  cut/2tt, 
starting  with  0.02  and  proceading  to  0.80  in  increments  of  0.02.  The 
third  card  producea  a  calculation  for  20  valuea  of  absciasa,  from  1.00 
to  1,95  in  atepa  of  0.05. 

The  program  printa  out  the  input  parametcra  for  reference  and  liats 
the  roaults  in  15  colurana .  The  firat  column  ia  the  abaciasa,  oit/Ztt. 

The  next  group  of  aix  columns  are  reaulta  for  the  fundamental  mode  only 
and  consist  of  Al,  Cl,  A2,  Al,  Cl,  and  A2.  The  next  aix  columna  are 
results  t*vr  nil  modea  up  to  11,  11.  Tha  printout  ia  in  the  order  A,  C, 
AFR,  B,  E,  MOMFR.  The  final  two  columns,  headed  X(z)  and  X(M)  indicate 
the  value  of  urt  at  which  the  negative  maximum  occurs  for  deflection  sr.d 
moment,  respectively . 

Until  the  abaciasa  exceeda  0.5G,  the  maximum  negative  responso 
during  forcing  is  zero.  The  program  doea  not  correctly  provide  for  thia 
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FREE  VIBRATION  MAXIMUM 

■C'MMU 


FIG.  A.  2  FLOW  CHARTS  FOR  PORTIONS  OF  ThE  MAIN  MULT.'PLATE  PROGRAM 

FOR  COMPUTATION  OF  DEFLECTIONS  (portion*  for  colculating  moment*  ore  similar) 
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□ 


SUBROUTINE  FREEZM 


!  ion  mu  »  ti 

I  wi  II  I 

("mwi  | 

SUBROUTINE  MODIFY 


A.3  SUBROUTINES  USED  IN  MULTIPLATE  PROGRAM 


S 


eventuality  sc  that  email  positive  or  negative  nWiA-ers  are  listed  for  C 
and  E  with  0  <  xt/2tt  <  0.50.  Such  values  should  be  disregarded. 

The  compile?  time  was  about  80  seconds  on  either  the  CDC3200  or  the 
B5500.  The  execution  time  per  abscissa  value  was  9  to  10  seconds  on  the 
3200  and  12  to  13  on  the  5500, 

The  results  of  the  program  are  depicted  in  Figs.  4,  5,  and  6,  The 
first  ia  the  maximum  deflection  or  moment  in  the  first  mode.  This  figure 
shows  the  fundamental  pattern  of  the  response  which  ia  only  slightly  modi¬ 
fied  by  the  participation  of  higher  modes.  The  first  maximum  to  occur 
during  forcing  is  celled  the  positive  maximum,  corresponding  to  the  posi¬ 
tive  pressure  of  the  N-wave.  The  curve  of  positive  maximum  is  monotoni- 
celly  increasing  with  sbsciasa  values,  approaching  20  at  infinity.  The 
negetive  maxima  occur  at  the  largest  negative  oscillation  during  the 
time  of  forcing.  Negative  maxima  may  occur  under  two  conditions;  et  the 
peak  of  an  oscillation  where  dw/dt  =  0  and  at  the  end  of  forcing  when 
t  =  t.  The  decreesing  portions  of  the  curve  pertain  to  the  firat  condi¬ 
tion,  increasing  portions  to  the  second.  The  free  maxima  show  the  lar¬ 
gest  values  end  moat  cleerly  show  the  dependence  of  the  mexima  on  the 
natural  frequencies.  For  sbaciaaa  values  of  1.43,  2,46,  and  3,47,  there 
is  no  motion  efter  the  time  of  forcing. 

The  multimodel  deflection  is  very  similar  to  the  first  mode  deflec¬ 
tion.  The  positive  maxima  show  e  2$  or  3$  increase  over  the  response  in 
the  fundamental  mode.  Negative  maxima  ere  very  similar  to  the  first  mode 
response  throughout  and  coincide  exactly  during  decreasing  portions.  The 
curves  for  fres  vibration  maxima  ahowa  slightly  broader  humps  and  sbouL 
a  4$  increese  for  the  first  hump,  less  for  later  humps.  The  minimum 
points  on  the  free  vibration  curve  are  not  zero,  showing  that  there  is 
always  some  residual  motion  following  forcing. 

In  Fig.  6  it  is  evident  that  the  moment  maxima  receive  a  significant 
contribution  from  the  higher  modes.  The  curve  of  positive  maximum  is 
augmented  up  to  20$.  In  spite  of  this  feirly  large  contribution,  the 
response  curve  is  smooth  except  for  a  slight  undulation  near  the  origin. 
The  average  increase  in  amplif icetion  factor  is  0.26,  The  curve  of  nega¬ 
tive  maximum  is  ?bout  5$  higher  than  thoae  from  the  first  mode.  The 
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curve  ia  bumpy,  showing  the  participation  of  higher  modea,  particularly 
the  aecond  symmetric  (1,3).  The  free  vibration  maxima  are  modified  moat 
algnlf icantlv  by  the  higher  modes.  The  contribution  of  the  1,3  mode  Is 
clearly  evidenced  by  the  triple  humps  at  each  main  hump  of  the  curve. 

The  increaae  over  the  first  mode  curve  varies  from  0  th  0.56  with  an 


average  value  of  0.28, 


PROGRAM 


MULTIPLATE 


C  PROGRAM  COMPUTES  CENTRAL  STRESS  ANO  DEFLECTION  OF  A  SQUARE,  PLATE  U'rCR 
C  OF  AN  N-WAVE.  FIRST  ELEVEN  MODES  'LINEAR  THEORY)  ARE  USED. 

COMMON  W  »D1 »X0 . J1 «K. OTAU  . ASIN*  . ACOSW 

DIMENSION  W (6 .6) .Z(4Q ) »R<6»6)»S(6.6) .MOM (40) «ASINW<6i6) » ACOSW < 6.6) 
HEAL  MO. MOM.  MOMFR . 1NCREM 

960  FORMAT (//8H  GROUP  . I2.6X.32HABSCISSA  =  OMEGA  TAU/2PI  FROM  . 

1  F6.J.UH  TO  »c6.3.10X.£HK=»I3»6X»4HEPS=.F6.3//> 

961  FORMAT  (11X.29HP  IRST  MODE  ONL  Y.19X.23HE  L  E  V  £  N 
1M  0  D  E  S) 

962  FORMAT  (9H  ABSC ISSA.2X.2HZ+.6X. 2HZ-.4X.6HZ  FREE.4X.2HM+.6X.2HM-. 

1  4X.6HM  FREE .  4X .  2HZ+ .  6X  .  2HZ- .  4X » 6HZ  FREE.4X.2HM-i-.6X.2HM-.4X. 

2  6HM  FREE. 4X.4HX(ZI »3X.4HX(MJ ) 

963  FORMAT  < 1X.F6.3* 12F8.4.2X. 2F7.3) 

980  FCRMAT(1X.13HSTATIC  DErL  S.F8.4.9X. 1SHSTATIC  MOMENT  =»F3.4) 

986  FORMAT  ( 13. I3.F6.3.F6.3.F6.3) 

987  FORMAT  €13: 

C  3EGIN  LOOP  TO  CALCULATE  ZO  AND  MO.  STATIC  DEFL  AND  MOMENT 
Z0=0 
M0=0 

DO  120  M=1.6 
FM=2*M-1 

R(M.M)=0.2S/FM**6 

S(M»M)=0.25/FM**4 

Z0-R(M»M)+Z0 

M0=S(M.M»*M0 

MM=M+1 

DO  11S  N=MM.6 
FN=2*N-1 

R(M,N)S(2.0/CFM*FN*(FM**2+FNA*2»**2) )*<-!. 0)<M  (M-N) 
S(M.N)=<1.0/(FM*FN*(FM**£+FN**2)  ))*<-1.0)**  <M-N) 

ZO=R(M.N)+ZO 
M0-SCMrN)+M0 
11S  CONTINUE 
120  CONTINUE 

WRITE  (6.980)ZC.M0 

C  CALCULATIONS  ARE  MADE  FOR  GROUPS  OF  VALUES  OF  ABSCISSA.  OMEGA  TAU/2PI 
C  ON  FIRST  DATA  CARD  ENTER  VALUE  FOR  THE  NUMBER  OF  DATA  GROUPS .JGROUP. 

C  ON  SECONO  AND  SUCCEEDING  DATA  CARDS  ENTER  JEND.  K.  FIRST.  INCREM. 

C  ANO  EPS  .  JEND  IS  THE  NUMBER  OF  INCREMENTS  OF  SIZE  INCPEM  IN  EACH 
C  GROUP.  K  IS  THE  TOTAL.  NUMBER  OF  STEPS  TAKEN  ON  BOTH  SIDES  OF  THE 

C  FIRST  GUESS. XO.  K  MUST  BE  000.  FIRST  IS  THE  FIRST  VALUE  OF 
C  ABSCISSA  IN  EACH  GROUP.  INCREM  DEFINES  THE  INCREMENTAL  CHANGE 
C  IN  ABSCISSA  SITHIN  EACH  GROUP.  EPS  IS  THE  ACCURACY  REQUIREMENT 
C  USED  IN  DECIDING  WHETHER  TO  ACCEPT  THE  CURRENT  MAX  OR  TO  RECYCLE 
READ  (5.987)  JGROUP 
DO  691  JUG=1. JGROUP 
READ  (5.986)  JEND» K.Fi  sT.INCREM.EPS 
FK=K 

Dsl.6/(FK-1.0) 

DO  690  111=1. JEND 
FIIUIII 

*eSCIS=FIRST*INCREM*(FIII-1.0> 


ACTION 

FORM  960 

FORM  961 

FORM  962 
FORM  962 


WRITE 
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OT  AU= ABSC I S*6 • 28  32 
00  iSO  M=l»6 
00  170  N=M»6 
FM=2*M-1 
FN=2*N-I 

W(MiN)=0TAU*(FM**2*FN*42)/2«  0 
ASINW(M*N)=SIN(W(M** ) ) 

ACOSW(M*N)=COS(W(M'N)) 

170  CONTINUE 
180  CONTINUE 

C  BEGIN  CALC  FOR  POSITIVE  OAF  FOR  DEFL  —  A  AND  A1 
Ol-D 

X0=2.0*Ai *l(0TAU/2.0) 

A1=1.0-CO.(XO  J-2.0*XO/OTAU+2.0*SIN(XC1)/OTAU 
C  CALC  9  VALUES  Ur  Z 
PENULTS 0 

250  CALL  FQRC74  (Z,R) 

CALL  ZMAX  17.) 

CALL  MOO IF 

CHANGES ABS ^ENULT-Z U 1 ) ) 

PENULT=Z (o 1  * 

IF  (CHANGE  2B0* 280.250 

200  A=Z(J1)/Z0 

C  BEGIN  CALC  FOR  POSITIVE  OAF  FOR  MOMENT *  B 
01=0 

X0=2 . 0* ATAN ! OTAU/2 . 0 ) 

PENULTsO 

310  CALL  FORCZM  (MOM.S) 

CALL  ZMAX  (MOM) 

CALL  MOOIFY 

CHANGEsABS( PENULT-MOM ( Jl) I 
P ENULT=M0M(J1> 

IF  (CHANGE-EPS)  340*340*310 
340  a  S  MOM(J1>/MO 
C  CALCULATE  OAF  NEG  FOR  Z 
UJJ  s  ABSC IS 
OJ JJsjjj 

XO  s  6.2832  •  OJJJ 

Cl  s  1.0  -  COS(XO)  -  ?.0  •  XO/OTAU  ♦  2.0  *SIN(XO)/OTAU 
IF  (JJJ-I)  402*402*404 
402  OI=0.5*OTAU/(FK-I .0) 

X0=0.75*OTAU 
GO  TO  410 

404  01= ( OTAU-6 .2832*0 JJJ ) / ( FK-3 • 0 ) 

IF  (01-01  406*406*408 
406  01=0 

40B  XO=OTAU/2«  0*3. I416*DJJJ-DI 

410  PENULTSO 

420  CALL  FORCZM  (Z*R) 

DO  430  1=1 • K 
ZlI)=-Z(I) 

430  CONTINUE 

CALL  ZMAX  (Z) 

CALL  MODIFY 

CHANGEsABS(PENULT-Z(  Jl) ) 


CALL  SUB 
CALL  SUB 
CALL  SUB 


CALL  SUB 
CALL  SUB 
CALL  SUB 


CALL  SUB 


CALl  SUB 
CALL  SUB 


55 


mn 


.  iiL  ii.ui-jj».ian.»"U"i «»» 


PENULi=Z<Jl) 

IF  (CHANGE-EP5)  455*455*420 
455  C=Z(J1)/Z0 

X0Z=XQ 

C  CALCULATE  OAF  neg  FOR  MOMENT 
IF(JJJ-l)  462* 462*464 
462  01=O.5*OTAU/tFK-1.0) 

XO=0*75*OTAU 
GO  TO  470 

464  0 is ( OTAU-6. 2B32*0 JJJ ) / (FK-3 . 0 ) 

IF  <0I-0>  466*466 * 46B 
466  01=0 

46B  XO=OTAU/2.1)+3.14I6*OJUJ”OI 

470  PENULTsO 

4B0  CALL  FORCZM  (MOM, 5) 

00  490  I=1*X 
MOM ( I >  s-MLM ( I ) 

490  CONTINUE 

CALL  ZMAX  (MOM) 

CALL  MOOIFY 

CHANGESAB5( PENULT-MOM i J1 ) ) 

PENULT=M0M(J1) 

IF  (CHANGE-EP5)  510*510*4B0 
510  E  =  M0M(JU/M0 
XOMsXO 

C  CALC  MAX  OEFL  DURING  FREE  VIBRATION 
01=0 

XQ=  0TAU/2.0  +  3.1416*(0JJJ+1.0) 

A 20=  -1.0  -COS(OTAU)  +2. 0*SIN(0TAU)/OTAU 
A2V=  SIN(OTAU)-  2.0/OTAU  *2.0*  C05(0TAU)/0TAU 
A2  =  A20  *C0S(X0-0TAU)  +  A2V  *SIN(X0  -OTAU) 
PENULT=0 

606  CALL  FREEZM  (Z*R> 

00  610  1=1 *K 
Z(II=ABS(ZU)» 

610  CONTINUE 

CALL  ZMAX  (Z) 

CALL  MOOIFY 

CHANGE=ABS (PENULT-Z ( J1 ) ) 

PENUL7=Z(J1) 

IF  (CKANGE-EP5)  620*620*606 
620  AFR=  Z(JD/ZO 

C  CALC  MAX  MOMENT  OURING  FREE  VIBRATION 
01=  0 

X0=  0TAU/2.0  ♦  3.1416* (OJJJ+1 .0) 

PENULTSO 

635  CALL  FREEZM  (MOM,S> 

00  640  1=1, X 
uryi  <  t  t  =  ABS  (  mom  i  I )  I 
640  CONTINUE 

CALL  ZMAX  (MOM) 

CALL  MOOIFY 

CHANGE=ABS (PENULT-MOM ( dl ) ) 

PENULT=M0M(J1) 

IF  (CHANGE-EPS)  650*650*635 


CALL  5UB 


CALL  5UB 
CALL  5UB 


CALL  SUB 


CALL  SUB 
CALL  5UB 


CALL  SUB 


CALL  5UB 
CALL  5UB 
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o50  MOMFRs  MOH(J1)/HO 

IF  (111*1)  660 * 660 * 661 

660  OATENO=JENO 
fcast=first*increm*(oateno-i.o> 

WRITE  (6*960)  JUG* FIRST* Ft. AST. It. EPS 
WRITE  (6*961) 

WRITE  (6*962) 

661  WRITE  (6*963)  ABSCIS* Al.Cl.A2* A1.C1.A2.A.C* AFR.B*E.MOMFR*XOZ*XOM 

690  CONTINUE 

691  CONTINUE 
ENO 

SUBROUTINE  FREEZM  (Y.G)  SUB 

COMMON  W.Ol.XO* J1*K*0TAU,ASINW*AC0SW 

OIMENSION  W(6*6)  *YU0)  *G(6*6)  .AS1NW(6*6)  *AC0SW(6*6> 

00  712  J=1*K 
01FFsJ-(K+l).'2 
XsXO  +01FF  *  01 
Y(J)  =  0 

7MX-0TAU)  712*704.704 
704  00  711  M=l*6 

00  7)0  N=M*6 
FM=2*M-1 
FN=2*N-1 

V=X*(FM**2+FN**2)/2.0 

Y(J>=  G(M.N)*( (2.0*ASINW(M*N)/W(M*N)-1.0-ACOSW(M*N) )*COS(V-W(M.N> ) 
1+(AS1NW(M.N)-2*0/W(M*N)+2,0*ACOSW(M*N)/W(M>N) )*S1N(V-W(M.N)))+Y(J) 

710  CONTINUE 

711  CONTINUE 

712  CONTINUE 
RETURN 
ENO 

SUBROUTINE  FORCZM  (Y.'O  SUB 

COMMON  W.Ol.XO* Jl.X.OTAU 
DIMENSION  W(6»6)*Y( 40) *6(6*6) 

00  B20  J=1*K 
UlFF=J-<K+l>/2 
XSXO  ♦  OIFF  •  01 
Y(J)  =  0 

1FIX-0TAU)  802*602*620 
802  IF ( X)  820*804.604 
804  00  619  M=l*6 

00  filfi  NSM.6 
FM=2*M-1 
FN=2*N-1 

V=X*(FM**2»FN**2)/2.0 

Y(J)=G(  I.N)*(1.0-COS(V)-2.0*V/W(M.N)*2.0*SIN(V>/W(M.N) >*YIJ> 
filfi  CONTINUE 
819  CONTINUE 
620  CONTINUE 
RETURN 
ENO 

SUBROUTINE  ZMAX  (Y)  SUB 
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COMMON  WiQl'XOtJIrKtOTAU 
OIMENSION  W<6*&»  «YUO>  .G<6.6> 
dl  =  I 
KKI=K-1 

00  870  J=I.KKI 

IF  <Y(J+I>-Y(JI>>  870.865. 86S 
665  JI  =J+I 
870  CONTINUE 
C  MAXIMUM  Y  IS  AT  Y<J1> 

RETURN 

END 

SUBROUTINE  MODIFY 
COMMON  H.01.X0.J1.K.0TAU 
DIMENSION  W(6»6> 

DIFF  =dl-(K+l>/2 

X0=X0+D1FF*D1 

FX-K 

0i=2.0*0I/<FK*I.0> 

RETURN 

END 
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APPENDIX  B 


RESPONSE  TO  INTERNAL  PRESSURE 


The  response  of  an  enclosing  structure  to  the  combination  of  boom 
loading  and  internal  pressure  will  he  considered  in  two  ways.  First,  it 
will  be  BBSumed  that  the  internal  pressure  ia  a  function  of  the  deflec¬ 
tion  of  the  element  considered.  Thia  corresponds  to  the  case  in  which 
the  moving  element  is  the  main  contributor  to  the  deflection  which  is 
ceusing  the  internel  pressure  rise.  For  the  second  celculatlon,  it  is 
assumed  that  the  internel  pressure  is  known  snd  is  not  modified  hy  motion 
of  the  element.  This  situetijn  occurs  if  the  motion  ef  the  element  is 
making  only  s  small  contribution  to  the  internal  prsseure. 

For  the  first  cese  the  relevant  equation  of  motion  is 


+  2 


a4w 

ftx^dy^* 


V  dsw 

g  bP 


q(t,x,y)  -  Gw  (B.l) 


where 

D  is  the  piste  stiffness 

w  is  the  deflection 

x,y  ere  rectangular  coordinates  in  the  plete 
V  is  the  weight/unit  eres 

g  is  the  acceleration  of  gravity 

t  is  t lme 

q  is  the  boom  loading 

G  is  e  coefficient  relating  deflection  end  internel  presses 

generated  by  ths  deflection.  Possible  units  ere  pounds/cubic 
inch. 

Following  the  derivetion  in  Chapter  5  of  Norris,  et  el.21,  let 

w  “  £  E  Dmn(t)  *mn  (B*2) 

end 

q(t,x,y)  =  f(t)  q(x,y)  qi 
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where 


D  (t) 

ion 

ia 

the 

^mn 

la 

the 

v(x>y) 

la 

the 

f(t) 

ia 

the 

q(x,y) 

is 

the 

Qi 

is 

the 

dynamic  ampiif icatlon  factor  for  each  mode, 

participation  factor  for  etch  mode 

ahape  factor  for  each  mode 

time  variution  of  applied  pressure 
spatial  variation  of  pressure 
peak  amplitude  of  pressure. 


ran 


A  uniform  pressure  distribution  over  the  plate  is  assumed,  so  that 
q(x,y)  s  1.0.  Then,  for  the  modal  deflection 

cf>mn(x,y)  =  sin  (mnx/a)  sin  (nTTy/a)  (B.3) 

the  participation  factor  is 

^mn  ~  TT8Dmn  ^ma+na)  a  for 
=  0  for  m  and/or 


m  and  n  odd  (B.4) 

n  even . 


The  participation  factcr  ia  identical  for  static  and  dynamic  loading. 
In  Norria,  et  al.ai,  it  is  shown  that 


00  00 


qta 


(x>y)  =  £  £  ♦, 


UJ 


v 

-  cp  (x,y) 

or  '  mn  '  f 


m=l  n=l 


mn  mn  g  Tmn 


(B.5) 


and  that 


r  +  — ' rr--;  +  t— r)  Cp  =  <jj*  —  cp 
'ax*  dxa3ya  3y  Yran  mn  g  Ymn 


(B.6) 


Now  Eq.  B ,2  is  substituted  into  Eq ,  B.l  taking  into  account  Eqa.  B.5 
and  B.6.  The  following  result  is  obtained: 


00  GO 


£  £  -  ♦  cp  +  <ua  D  +  — £  D  -  UJa  f(t)]  =  0  (b.7) 

t-*  g  Trnn  Trnn  *■  d  t a  ran  mn  v  mn  mn  '  ' J 


mal  n=l 


mn  mn  y  mn  mn 
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In  order  to  satisfy  this  equation,  eech  of  the  terns  of  the  doubly  infinite 
series  must  be  zero.  The  solution  for  eny  term  is 


mn 


t 

7/7*0.  I  t(t')  sin  l/ul^  +  d»/V  (t-f 
mn  — *  hd 

v 


)3  dt*  (B.8) 


Compere  Eq.  B.8  witb  the  solution  for  external  loading  only,  wbich 
is  t 

(B.8) 


D  =  ci) 
mn  mn 


f(t')  sin  [uL_(t-t')]  dt' 
mn 


The  effect  of  internal  pressure  Is  to  increese  the  frequency  by  the  retio 

/u,2mn+ 


Rjj(n>,n)  = 


IB 


( B ,  10  ) 


mn 


end  to  decreese  the  dynamic  emplif icetion  factor  by  the  seme  ratio.  For 
the  further  discussion  only  the  first  mode,  cb1x,  will  be  considered. 

For  this  cese  Fig.  4  cen  be  used  to  determine  D1X  provided  the  frequency 
end  DX1  ere  modified  by  the  retio  in  Eq .  B.10. 

For  the  second  cese  the  intemel  pressure  is  assumed  to  be  sinu¬ 
soidal,  given  by  the  equation 

q  =  q2  sin  2TTt/r  during  forcing 


=  0 


after  forcing 


T  is  the  duration  of  forcing,  end  coincides  with  the  boom  duration 
q 3  is  peek  applied  pressure. 

For  e  combination  of  boom  and  internal  pressure,  the  loading  is 


q  =  cii(  l-2t/T )  -  qa  sin  2nt/T  (t  *  t) 
=  0  (t  >  t) 


(b .11) 
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?  solution  was  obtained  for  the  first  mode  only  using  the  methods  of 
Appendix  A.  The  deflection  during  forcing  is 


Z  2 

—  =  qt  [l-cos  (tut )  -  2t/T  +  —  sin  (wt)] 

(B. 12) 

-  l^TT3/^Ta  C+  8in  (^/t)  -  2tt/wt)  sin(ujt)3 

where 

Z1  is  the  static  deflection  under  qj  =  1.0 

With  internal  pressure  slone  (qj  =  0)  the  maxima  occur  for 


COS  (2TTt/T)  =  COS  (rjut) 
The  relevant  maxima  occur  at  the  times 


t  = 


2nrr 


2tt 

“  +  T 


(B.13) 


The  time  of  positive  maximum  is  given  by  n  ^  1.  The  time  of  negative 
maximum  is  given  by  the  largest  even  vslue  of  n  such  thst  t  £  t. 

During  free  vibration  the  deflection  is 


Z  I  2 

=  qX  [-1.0  -  cos  (u/r)  +  —  sin  (pjt)]  cos  [u)(t  — r) D  +  [“sin  (out) 


Z1 


”  +  ”  coa  (u^)3  sin  U(t-T)]}  -  qa  [1_^a(ly4TTa]  (  sin  (wt)  cos  [u>(t-T)j 


+  [-1.0  +  cos  (u>t)]  sin  [iu( 


f 

t-T)]^ 


(B.14) 


For  q}  =  0,  the  condition  for  a  maximum  is  that 


COS  ['ju(  t-T )  3  =  cos  (ujt) 


That  is,  for 


Pit  n  1  qT 

2tt  =  2  2  2n 


(B.15) 


62 


Two  programs  were  written  to  carry  out  the  required  computstione . 
The  one  called  Boom-and -Internal  makee  calculations  of  maxima  when  both 
qt  and  q2  are  not  zero.  No  flow  chart  ia  provided  for  it  ss  the  program 
ie  a  aimplif icstion  of  the  Multipiste  program  in  Appendix  A.  The  input 
parameters  are  the  same  as  for  Multiplate  with  the  addition  of  Qj  and  Q2 
which  are  the  coefficiente  of  the  boom  and  internal  preeaure  terms.  The 
maxima  of  Eqs.  B.12  and  B.14  were  determined  by  trial  using  the  same 
eyatem  as  in  Multipiste. 

When  qt  =  l=  0,  the  terme  of  the  maxima  can  be  determined  analy¬ 
tically,  so  «  special  program  called  Sinusoidal  was  written  to  calculate 
this  case.  There  are  vto  data  cards  needed  for  this  program  aa  the  input 
is  part  of  the  program  instructions. 

Both  programs  were  written  in  Fortran  IV  and  run  on  a  Burroughs 
B5500  computer.  Compilation  times  were  about  50  seconds.  Sinusoidal 
executed  at  a  rate  of  0.12  seconds  per  sbeclsss  value,  Boom -end -Internal 
took  0.15  seconds  for  each. 

The  resulta  of  the  Sinusoidal  program  are  in  Fig.  7.  This  shows 
some  of  the  same  features  ss  the  response  to  booms  but  the  first  peaks 
are  higher  and  later  peaks  are  lower. 

The  resulte  of  Boom-and -Internal  for  Q2  =  0.25  and  0.50  are  ebown 
in  Figs.  8  and  9.  Evidently  the  main  effect  of  Internal  pressure  is  to 
decrease  the  firet  hump  in  the  response  curve.  Otherwise  the  curves  are 
essentially  the  esme  ss  those  in  Fig.  4  for  boom  loading  only. 
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PROGRAM 


BOOM  AND  INTERNAL 


C  PROGRAM  COMPUTES  MOTION  OF  A  1  DOF  SYSTEM  TO  A  BOOM  LOADING  AND  AN 
C  INTERNAL  PRESSURE  IN  THE  FORM  05  A  SINE  WAVE 
COMMON  D1 *  XO * J1 * K • OT AU ♦ Z * A8SC 15 • »1 * 02 • EPS 
DIMENSION  Z(40J 

901  FORMAT  <9H  ABSC ISSA  *  4X*  4HDP0S * 4'X * 4HDNEC *  3Y  * 5HDFREE  *  4X  * 4HTP0S *  4X  * 

1  4HTNEG*  3X * 5HTFREE ) 

9S0  FORMAT l S2H  DAF  FOR  1D0F  SYSTEM  UMOER  SINUSOID  AND  800M  LOADING) 
955  FORMAT  </10X*2HKs* 12*4X*4HEPS=*F6.3*4X*3HQl=*F6.3*4X*3HQ2=*F6.3/> 

963  FORMAT  ( 1X*F6.2*6F8.3) 

986  FORMAT  ( 13* I3*SF6.3 ) 

987  FORMAT  (13) 

WRITE  (6 #950) 

READ  (5# 987)  JGROUP 
DO  691  JUGs 1 ► JGROUP 

READ  (5r986)  JEND#K»F1RST» DELTA#EPS»Q1 #02 
WRITE  (6#9SS>  K#EPS#01#02 
WRITE  (6*901 ) 
fk=k 

D=l.b/(FK-1.0) 

DO  690  111=1* JEND 
F111=III 

ABSCIS  =  FIRST +OELTA  ♦(F1II-1.0) 

OTAU  =  ABSCIS*6.2832 
C**a**C ALCULATE  ZPOS 
D1=D 

X0=2.0*ATAN(OTAU/2.O) 

PENULT  =  0 
250  CALL  FORCEZ 
CALL  ZMAX 

CHANGE  s  ABS ( PENULT-Z ( J1 ) ) 

PENULT  =  Z(J1) 

IF  (CHANGE-EPS)  280*280*250 
280  ZPOSsZ(Jl) 

XPOS  =  XO/OTAU 
C*****CALCULATE  zneg 
Jjj  =  ABSCIS 
DJJJsJJJ 

XO  =  6. 2832*0 JJJ 
IF  (JJJ-1)  402*402*404 
402  Dl=0.S*0TAU/(FK-1.0) 

X0=0.7S*0TAU 
GO  TO  418 

404  01= ( OTAU-6 . 2832*DJ JJ  > / (FK-3 . 0 ) 

IF  (Dl-D)  406*406*408 
406  Dl=D 

408  X0=0TAU/2.0+3.1416*DJJJ-Dl 

91s  FZNJL7  —  C 
420  CALL  FORCEZ 
DO  430  I=i*K 
Z ( I )  =  -Z(I) 

430  CONTINUE 
CALL  ZMAX 
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CHANGE  =  ABSIPEN'jlWUD) 

PENULT=Z(J1> 

IF  (CHANGE-EPS)  455.455.420 
455  ZNEG  =  Z(J1> 

XNEG  =  XO/OTAU 

c*****calculate  zfree 
01=0 

XO  =  OTAU/2.0+3.I4I6*(OJJJ+l*O> 

PENULT  =  0 
606  CALL  FREEZ 

00  610  I=I»K 
Z(I)  =  ABSlZlin 
610  CONTINUE 
CALL  ZMAX 

CHANGE  =  ABSIPENULT-Z(JI) > 

PENULT=Z( JI> 

IF  (CHANGE-EPS)  620.620.606 
620  ZFREE  =  Z(UI> 

XFREE  =  XO/OTAU 

WRITE  (6.963)  ABSCIS.ZPOS.ZNEG. ZFREE. XPOS.XNEB.XFREE 

690  CONTINUE 

691  CONTINUE 
END 

C 

C******SUBROUTINE  FORCEZ 
SUBROUTINE  FORCEZ 

COMMON  01. XO.JI.K.OTAU. Z.L3SCIS.QI.Q2.EPS 
DIMENSION  Z (40 > 

00  820  J=I»K 
OIFF  =  J-(K+I)/2 
X  =  XO  ♦  OIFF  *  01 
Z(J>  =  0 

IF(X-OTAU)  B02.802.B20 

802  IF(X)  B20.B03.B03 

803  AA=ABS(ABSCIS-I.O) 

ZI  =  1.0-C0S(X)-2.0*X/OTAU*2.0/OTAU*SIN(X> 

IF  (AA-EPS)  8I0.B10.B04 

B04  Z2=  -ABSCIS**2/(1.0-ABSCIS**2)*t-SIN(X/ABSClS>  ♦  SIN (X) /ABSCIS) 

GO  TO  BI5 

810  Z2=  -0.5*(SIN(X)-X*C0S(X>  > 

BI5  Z(J)  =  QI*Z 1  ♦  Q£*Z2 

820  CONTINUE 
RETURN 
END 

. . ^SUBROUTINE  FREEZ 

SUBROUTINE  FREEZ 

COMMON  01.  XG.J1.K.0TAU. Z.ABSCIS.GI.Q2.EPS 
DIMENSION  Z (40 ) 

OO  712  J=1.K 
OIFF  =  J-(K*l>/2 
X  =  XO  +  OIFF*OI 
Z(Ji  =  0 

IF (X-OTAU)  712.704.704 
704  AA=ABS( ABSCI5-1 .0) 

Zl  =  ( -I « 0~C9S ( OT  AU ) +2 • 0*SIN ( OT  AU) /O  T  AU } *COS ( X  -OTAU) 
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1  ♦(SIN(OTAU)-2.n/OTAU*2.0*COS(OTAU)/OTAU)*SIN(X  -GTAU) 

IF  (AA-EPS)  710.710*706 

706  Z2=-(ABSCIS/(1.0-ABSCIS**2>)*(SIN(0TAU)*C0S(X  -OTAU) +(-1.0+ 

1  COS (OTAU) ) *SZN<X  -OTAU)) 

SO  TO  711 

710  Z2=  3.1416«C0S ( X-OTAU) 

711  Z(J)  =  Q1*Z1  ♦  Q2*Z2 

712  CONTINUE 
RETURN 
END 

^♦••♦♦SUBROUTINE  ZMAX 
SUBROUTINE  ZHAX 

COMMON  Dl.XO.dl.K. OTAU. Z.ABSCIS.Q1.G2.EPS 
DIMENSION  Z (40 ) 

Jl=l 

KKlsK-1 

UO  870  d=l.KKl 

IF  <Z(d+l)-Z<dl) )  B70.86S.B65 
B65  dl  =  dn 
870  CONTINUE 

01FF  S  dl-(KM>/2 
XO  s  XO  ♦  D1FF  •  Dl 
FK  —  K 

Dl  =  2.0*01/ ( FK+1 .0) 

RETURN 

END 

?DATA  C07BSEAMAN/TAPE5 


4 


25 

21 

0.5 

0.02 

0.005 

1.0 

0.0 

20 

21 

1.6 

0.02 

0.005 

1.0 

0.0 

15 

21 

2.7 

0.02 

0.005 

1.0 

0.0 

15 

21 

3.7 

0.02 

0.005 

1.0 

0.0 

7EN0  OF  DECK 
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PROGRAM  SINUSOID 


C  SINUSOID  CALCULATES  THE  PEAK  POSITIVE *  NEGATIVE*  AND  FREE  RESPONSE 
C  OF  A  1  DOF  SYSTEM  TO  A  LOADING  IN  THE  FORM  OF  A  SINE  WAVE  EXTENDING 
C  FROM  ZERO  TO  2  PI  i  WITH  NO  LOADING  THEREAFTER 

900  FORMAT  (1X*F6.2*2X *6F8. 3) 

901  FORMAT  ( 9H  ABSCISSA*4X*4HDP0S*4X*4HDNEG*3X*SHDFREE*4X*4HTP0S.4X* 
1  4HTNEG*  3X * 5HTFREE ) 

902  FORMAT  (35H  DAF  FOR  IDOF  SYSTEM  UNDER  SINUSOID) 

WRITE  (6*902) 

WRITE  (6*901) 

DO  100  11=1*200 
IF  (11-50)  10*5*10 
5  0P0S=1.5708 

DNEG=-3.1416 
DFREE=3« 1416 
TPOS=0.5 
TNEG=1.0 
TFREE=1.0 
ABSC1S  =  1.0 
GO  TO  90 
10  BIi=II 

ABSC1S=BII*0.02 

DPOS  =  ABSCIS/ ( I • 0-ABSCIS) *S IN (6. 2B32* ABSCIS/ ( 1 . 0+ ABSCIS) > 

TPOS  =  1 .0/ ( 1 • 0+ ABSCIS) 

M  s  0.5  ♦  0.5*ABSC1S 
FN  =  2*M 

UNEG  =  ABSCIS/ ( 1 .0-ABSCIS) *S IN (6 « 2832*FN* ABSCIS/ ( 1 • 0+ABSCIS) ) 
0NEG2=  ABSCIS/ (1.0-ABSCIS**2)*SIN(6.2')32*ABSCIS) 

TNEG  =  FN/ ( 1 . 0+ABSC IS) 

DIFF  =  -DNEG  ♦  DNEG2 
IF  (DIFF)  15*20*20 
15  DNEG=DNE62 

TNEG=1 • 0 

N  =  ABSCIS  ♦  1.0 

20  DFREE  =  2.0*ABSCIS/(1.0-ABSCIS**2)*SIN(3.I416*ABSCIS)»(-1.0)**N 

M  =  ABSCIS  ♦  1.5 
FM  =  2*  M-l 

TFREE  =  0.5  ♦  FM/4. 0/ABSCIS 

90  WRITE  (6*900)  ABSC IS* DPOS *DNEG* DFREE* TPOS*  TNEG* TFREE 

100  CONTINUE 

END 
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APPENDIX  C 


NONLINEAR  LOAD -DEFLECTION  RELATIONS 

The  purpose  of  this  appendix  is  to  determine  the  form  of  the  lerge- 
def lection  relations  between  load  and  deflection,  end  between  etress  and 
deflection  for  squere  pistes.  These  reletlone  will  be  used  In  the  next 
eppendlx  ss  e  bssls  for  reducing  evellable  experimental  date. 

In  the  litereture  the  lerge  deflection  of  square  pletes  hae  been 
determined  only  for  some  specific  values  of  Poisson's  ratio  end  certein 
boundery  conditions;  the  sveileble  solutions  do  not  precisely  fit  the 
present  requirements.  An  approximate  rsther  than  rigorous  solution  seemed 
adequate,  snd  It  was  also  desirable  to  have  the  reeults  in  analytical 
rather  than  grsphlcsl  form.  A  modlficetion  of  the  method  recommended  by 
Ftippl  (1924)  wss  used  for  the  calculations . 

The  method  of  Ftippl  Is  to  seek  e  solution  as  tbe  sum  of  two  parts: 
one  for  membrane  snd  one  for  bending  behavior.  Thus  tbe  applied  loading 
csn  be  represented  es 


«  -  9t,  +  «m 


(C.l) 


where 

qs  * 

Q  =  ,  s  nondime?'  lonsl  loading  peremeter  end  subscripts  refer 

to  bending  and  membrene. 

q  =  pressure 

s  =  length  of  the  eide  of  the  square  plete 

In  F&ppl’s  approach,  the  relations  between  the  Q's  end  the  deflec¬ 
tion  sre  ieken  from  the  solutions  of  the  small  deflection  bending  and  of 
the  membrane  problems.  These  solutions  have  the  form 


%  s  Ao?  * 

Qm  -  B0?» 

so  that  Eq.  C.l  becomes  pR£CED!KS  PAGE  BLANK 


(C.2) 
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Q 


(C.3) 


=  A0?  +  B0£3 


where 

5  is  the  nondimensions  1  deflection  w0/h 
w 0  is  the  central  deflection 
A0  and  B0  are  functions  of  Poisson's  ratio. 

When  thia  cubic  equation  is  solved,  the  membrane  and  bending  stresses 
can  be  determined  from 


Sb  =  Col 
S„  -  D,| 


(C.4) 


where 


a  =  stress 

C0  snd  D0  sre  functions  of  Poisson’s  ratio,  v. 

The  natural  boundary  conditions  for  this  solution  are  the  simple  support 
for  bending  and  the  immovable  edge  for  membrane  action.  However,  in  our 
case  the  edge  is  movable.  Therefore,  the  sought  solution  was  assumed 
to  have  the  form  of  Eqs.  C.3  and  C.4  but  with  unknown  coefficients  Bt 
and  D1.  Further,  it  wss  assumed  thst  B  and  Bx  would  vary  in  the  same 
way  ss  a  function  of  v  so  that 


hi 

(V.) 

Bi 

(C.5) 


snd  similarly  for  D  and  Di. 

The  steps  tsken  in  the  solution  were  the  following: 

1.  Find  the  relations  between  deflections  and  pressures  for  the 
separate  esses  of  raembrsne  snd  bending,  that  is,  find  A0(u) 
and  B0(uj  in  Eq.  C.2. 

2.  Use  the  theoreticsl  solution  of  Levy18  to  determine  the  coef¬ 
ficients  Bj  snd  Dx  (for  the  esse  u  =  0.316), 

3.  Determine  the  coefficients  Bx  and  Di  for  the  case  y  =  0,23 
UBing  Fq.  C.S. 
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The  values  of  A0  and  C0  sre  readily  determined  from  Chapter  5  of 
Timoshenko  (1959),  The  results  are  given  In  Table  C.l. 

TABLE  C,1 

COEFFICIENTS  FOR  LARGE  DEFLECTION  OF  SQUARE  PLATES 


V 

a 

Bo 

Co 

Do 

Bi 

Di 

0.23 

21.7 

28.65 

5.91 

2.134 

12.8 

1 

0.316 

22.8 

31.3 

6.635 

2.582 

14.0 

The  coefficients  B  snd  D  were  determined  from  the  snslyals  on  p.  419 
of  Timoshenko  (1959)  for  the  roembrsne  problem.  The  following  expressions 
were  taken  for  the  displacements 


w 

u 

V 


wo 

cos 

T7X 

s 

cos 

31 

s 

c 

sin 

2ttx 

a 

cos 

31 

a 

c 

a  in 

2ny 

coa 

nx 

(C.6) 


where  u,  v  are  displacements  in  the  plsne  of  the  piste  in  the  x  snd  y 
directions,  respectively,  and  the  origin  is  at  the  center  of  the  plate. 
The  strain  energy  of  the  piste  is  then 


V 


Eh  /  5tt4  4 
2  ( 1- v  2 )  \  Mi1  W° 


5-3v 

3a 


tt2c»03  + 


t2?  "2  ♦  f 


The  principle  ol  virtual  displacements  wss  used  then  as  outlined  in 
Timoshenko  to  find  c  snd  w0  ss  functions  of  the  applied  pressure, 
values  of  B  snd  D  in  Table  1  sre  the  calculatod  results. 


(C.7) 


The 


The  theoreticsl  values  of  Levy18  were  used  to  find  Bt  snd  Dx.  Hia 
results  were  plotted  e.a  membrane  and  bending  stresses  versus  Bpplled 
pressure  for  the  caf.e  v  =  0.316,  Two  relations  were  uaed  for  the  deter¬ 
mination.  The  first  was  the  requirement  that  the  deflection  in  Eq.  C.4 
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be  unique ,  that  is,  that 


(C.8) 


Levy's  results  were  used  for  S.  and  S  and  a  value  of  Dx  =  1.40  waa  found 

b  m 

(aee  Fig.  C.l).  Then  Bqa.  C.3  and  C.4  were  combined  aa 


q 


+  Bx 


/Sa,s/a 


{  C .  9  ) 


Using  Levy's  values  for  Q,  and  S^,  a  solution  waa  obtained  for  B1# 

The  determination  of  Bj  ia  shown  in  Fig.  C.2  from  which  it  ia  evident 
that  Eq.  C.9  ia  a  good  representation  of  the  relation.  Following  this 
evaluation  of  Bt  and  Dx  for  u  =  0.316,  Bj  and  Dt  for  v  0.23  were  cal¬ 
culated  using  Eq.  C.5. 

The  following  equations  are  then  available  for  solving  large  deflec¬ 
tion  problems  in  aquare  plates  with  aimply  supported  boundaries. 


Q  =  A05  +  Bj 

sb  =  C0§  (C.10) 


In  the  next  appendix  experimental  valuea  are  used  to  re-evaluate  the 
coefficients  A0  through  D1  in  these  equations. 
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FIG.  C.2  GRAPH  FOR  DETERMINATION  OF  B( 
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APPENDIX  D 


EXPERIMENTAL  DATA  ON  WINDOW  RESPONSE 

Load-daf lactlon  and  atrass-daf lactlon  ralatlor.a  can  ba  obtained 
from  raaulta  of  aavaral  experimental  studlaa.  Thla  atudy  la  basad  on 
tha  results  of  Bowlaa  and  Sugarman  (1952),  Orr  (1957),  and  Praynlk  (1963). 
Tha  data  ara  fitted  to  aquatlona  In  tha  form  of  thosa  in  Appendix  C. 

Load-Daf lactlon  Relations 

A  load-daf lactlon  aquation  Ilka  Eq.  C.3  will  be  datarmlnad  flrat. 

Q  a  A0?  +  B0f 3 


whara 

Q 

= 

nondl mens  Iona  1  load  =  §  (§•) 4 

E  n 

1 

= 

nondimans Iona 1  daflactlon  =  w 

o/h 

q 

= 

applied  pressure 

E 

= 

Young's  modulus 

a 

dlstanca  between  supports  of 

square  pana 

h 

= 

thlcknasa 

= 

center  deflection 

Tha  value  of  A0  will  be  takan  as  21.7,  tha  thaoratlcal  valua  obtained 
from  small-daflactlon  plata  theory  for  Polason’s  ratio  aqual  to  0.23. 

Tha  data  will  ba  ralatad  to  this  aquation  to  datarmlna  tha  valua  for 
B0,  and  to  Indicate  how  wall  tha  data  follow  tha  analytical  aquation. 

Wa  will  usa  tha  data  on  ultimata  deflection  (deflection  juat  pracadlng 
fractura)  tr  i  Orr  (1957)  and  from  Bow las  and  Sugarman  (1952).  Tha 
raaultlng  equation  should  be  aspaclally  applicable  at  tha  utllmate 
strength  of  panala,  but  will  not  necessarily  describe  the  load-daf lactlon 
relation  at  low  pressures .  Tha  data  used  ara  In  Tables  D.l  and  D.2. 

0  r's  ra.sults  represent  Individual  panela  wharaas  Bowlas  and  Sugarman 
tested  30  or  40  panels  of  aach  slza  and  reported  mean  valuas  and  stan¬ 
dard  davlstlons  from  the  maan.  Tha  raaults  are  traatad  hara  as  though 
thay  wara  from  square  panala  although  Orr'a  panels  wara  rectangular, 
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TABLE  D.l 


ULTIMATE  DEFLECTION  AND  BURSTING  PRESSURES 
(Dats  from  Orr) 


h 

Thickness 

_C*pA 

s  '* 

ares 

JinA) 

a/h 

q 

Msx .  Pres . 

_ Cep*} _ 

Q* 

Nondiin . 
Pressure 

wo 

Max,  Defl. 
(in.) 

e  -  to. 
s  h 

PLATE 

0.2373 

6720 

346 

52.25 

516 

1.200 

5.05 

0.240 

6720 

342 

51,87 

491 

1.189 

4.95 

0.303 

6720 

271 

80.65 

300 

1,200 

3.96 

0.301 

6720 

273 

56.17 

214 

1.000 

3.32 

0.2344 

8360 

389 

39.26 

631 

1.300 

5.54 

0.2453 

8360 

373 

36.02 

,84 

1.200 

4.89 

0,3045 

8360 

300 

54.09 

305 

1.200 

3.94 

0,305 

8360 

299 

54.02 

303 

1.200 

3.93 

0.242 

9840 

410 

32.52 

638 

1.400 

5.78 

0.239 

9840 

416 

23.58 

485 

1,200 

5.02 

0.303 

9840 

327 

44,56 

355 

1.311 

4.33 

0.304 

9840 

327 

44.00 

343 

1,300 

4.28 

0.369 

9840 

269 

56.13 

204 

1.200 

3.25 

0.372 

9840 

267 

57.85 

203 

1.200 

3.23 

0.114 

8640 

816 

16.72 

5120 

1.400 

12.3 

SOLEX 

D 

8840 

387 

54.22 

740 

1.400 

|  5.65 

ns 

8640 

257 

73.20 

197 

1.200 

3.22 

0 .355 

8640 

2a  1 

61.36 

148 

1.000 

2.61 

0.255 

11530 

421 

38.56 

847 

1 

5.88 

0.248 

11530 

433 

41.52 

1015 

1 

6.45 

*  Q  »  t  E  was  tsken  as  107  psi. 

CD 


ULTIMATE  DEFLECTION  AND  BURSTING  PRESSURES 


with  sspect  ratios  between  0.6; 1.0  and  1.0; 1.0.  The  dsta  points  sre 
plotted  in  Fig.  D.l.  The  trend  line  is  given  by  the  equation 

Q  =  21.75  +  2.805s  (D.l) 

which  sppears  to  fit  the  date  quite  well.  For  comparison  the  theoretical 
results  of  Levy  (1942)  for  Poisson's  rstio  of  0.316  are  also  shown. 

A  check  can  be  made  to  determine  how  well  Eq.  D.l  predicts  deflec¬ 
tions  st  low  loads,  using  the  data  of  Bowles  and  Sugsrman  on  four  panels 
(listed  in  Table  D.3).  These  four  sets  of  load-deflection  points  are 
shown  in  Fig.  D.2  with  the  theoreticsl  curve  of  Eq.  D.l.  Eq.  D.l  will  be 
used  in  subsequent  work,  with  the  understanding  that  its  accuracy  is 
limited  at  small  deflections.  It  may  be  noted  that  the  data  from  the 
tests  of  the  l/8-inch  piste  are  far  from  the  other  points.  This  tendency 
is  noted  in  later  graphs  also.  The  discrepancy  may  he  caused  hy  using 
the  mean  thickness  of  0,122  inch  in  cslculsting  Q,  and  S  or  there  may 
have  heen  some  experimental  error. 


FIG.  D.l  RELATION  BETWEEN  APPLIED  PRESSURE  AND  DEFLECTION  AT  RUPTURE 
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TABLE  D.3 


8TRE88  DATA 

(Date  from  Bowles  and  Sugerman) 


h* 

Thickness 

(in.) 

q 

Pressure 

(psl) 

Q** 

Nondlm. 

Pressure 

S  =  Wo/h 
Nondlm. 
Deflection 

°b 

Bending 

Stress 

(psl) 

am 

Membrene 

Stress 

(psl) 

sb** 

Nondlm. 

Bending 

Stress 

S  ** 
NonSlm. 
Membrene 
Stress 

0.122 

0.05 

57.8 

1.61 

595 

215 

6.40 

2.315 

0.10 

115.5 

2.40 

795 

415 

8.55 

4.465 

0.15 

173  .2 

2.98 

900 

580 

9.68 

6.24 

0.20 

231.0 

3.43 

950 

750 

10.23 

8.07 

0.25 

289.0 

3.78 

0.195 

0.05 

8.9 

0.59 

577 

92 

2.43 

.39 

0.10 

17.8 

0.92 

950 

200 

3.99 

.84 

0.15 

26.6 

1.19 

1222 

297 

5.14 

1.25 

0.20 

35  .5 

1.41 

1430 

400 

6.01 

1.68 

0.25 

44.4 

1.62 

1595 

525 

6,71 

2.21 

0.30 

53.3 

1.81 

1725 

645 

7.25 

2.71 

0.35 

62.2 

1.96 

1827 

752 

7.68 

3.16 

0.25 

0.1 

6.55 

0.42 

670 

40 

1,715 

0.10 

0.2 

13.10 

0.70 

1240 

160 

3.17 

0.41 

0.3 

19.65 

0.93 

1705 

295 

4,37 

0.75 

0.4 

26.20 

1.13 

2075 

435 

5.31 

1.11 

0.5 

32.75 

1.32 

2367 

563 

6.06 

1.44 

0.6 

39.30 

1,48 

2615 

685 

6.69 

1.75 

0.7 

45  .85 

1.63 

2825 

815 

7.23 

2.09 

0.8 

52.40 

1.76 

2990 

950 

7.65 

2.43 

0.373 

0.2 

2.65 

0.25 

595 

45 

0,684 

0.052 

0.4 

5.42 

0.41 

1170 

100 

1.345 

0.115 

0.6 

7.94 

0.57 

1730 

180 

1.99 

0.207 

0.8 

10.58 

0.70 

2270 

270 

2.61 

0.311 

1.0 

13,23 

0.81 

2740 

380 

3.15 

0.437 

1.2 

15.87 

0.92 

3170 

520 

3.64 

0.598 

1.4 

18.52 

1.03 

1.6 

21.15 

1.12 

|  1.8 

23  .80 

1.21 

2.0 

26  .47 

1.30 

2.2 

29.10 

1.38 

*  The  thickness  of  the  specimen  tested  was  not  known,  so  the  mesn 
thickness  for  the  group  of  specimens  was  used. 

**  Nondimensionel  pressures  snd  stresses  wore  celculeted  using  E  -  107  psl. 
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NONDIMENSIONAL  PRESSURE,  0 

t*-so*s-i« 

FIG.  D.2  EXPERIMENTAL  AND  THEORETICAL  LOAD-DEFLECTION  CURVES 
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Ultimate  Deflection  end  Ultimate  Strength 


The  deflection  which  e  penel  cen  undergo  before  falling  Is  closely 
releted  to  the  thickneaa  of  the  panel.  The  deflection  data  of  Tebles 
D.l  end  D.2  ere  plotted  In  Figs.  D.3  end  D.4  to  show  thla  reletlon.  It 
lfl  cleer  from  the  figures  thet  there  la  e  marked  difference  between  the 
reaulta  of  Orr  (1957)  and  of  Bowles  and  Sugarman  (1952),  end  there  is 
e  noticeeble  difference  between  tbe  plete  and  abeet  gless  reaulta  In  the 
lattor  reference. 

To  predict  ultimate  strength  of  gleaa  penes  It  la  useful  to  plot 
nondlmenaionel  maximum  preaaure  agelnat  the  thickneaa  retlo,  e/h,  es  In 
Fig.  D.5.  Ageln  there  is  e  distinct  difference  between  the  results  of 
Botflea  and  Sugerman  (1952)  end  of  Orr  (1957).  The  pressures  taken  by 
the  panes  of  the  former  ere  about  2.5  timea  thoae  In  the  letter.  This 
difference  Is  only  partly  eccounted  for  by  the  feet  that  Orr's  penea  were 
rectanguler  and  tested  8t  e  slower  loading  rete. 

To  Indicate  the  probable  effect  of  testing  rate  on  ultimate  strength, 
we  can  exemlne  the  equetlon  given  by  Frownfelter  (1959) 


o 

max 


+  0.5  log10 


da 

dt 


where 

da/dt  la  the  rete  of  eppllcatlon  of  atresa  in  pai/min 
a0  la  e  reference  atreas 

With  a  reference  atreas  of  2000  psl,  thla  equetlon  la  plotted  in  Fig.  D.6. 
In  addition,  the  recommended  strength  veluea  from  Pittsburgh  Plete  Gleas 
(1964)  ere  shown  and  straight  lines  were  drewn  on  the  figure  to  lndlcete 
the  trend  of  the  points,  Comparison  of  the  testing  retea  of  Bowles  and 
Sugarman  (1952)  with  thoae  of  Orr  (1957)  indicates  that  Orr'a  strength 
velues  should  be  ebout  75‘jfc  of  those  from  Bowles  end  Sugarman. 


Streaa-Def lection  Hatloa 

Equations  relating  maximum  atreaa  end  centrel  deflection  were  deter¬ 
mined  ualng  tl.2  date  of  Bowlea  end  Sugerman  (1952).  The  equetlona  have 
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I  I  I  I  I  I  I 

BOWLES  AND  SUGARMAN  DATA 


900 

•005-17 


FIG.  D.3  DEFLECTION  AT  RUPTURE  AS  A  FUNCTION  OF  THICKNESS  RATIO,  o/h 


RATIO,  c/h 

TB-C0«5-I9 


FIG.  D.5  FAILURE  PRESSURE  ON  PLATES  AS  A  FUNCTION 
OF  THE  THICKNESS  RATIO,  aA 


BREAKING  STRESS,  p* 
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the  same  form  as  those  of  Appendix  c,  that  Is 


Sb  =  Co? 
Sm  -  Do?’ 


where 


a  =  stress 
C0 ,  D0  =  constants 

b,m  =  subscripts  denoting  bending  or  membrane. 

The  values  used  ere  listed  In  Table  D.3  for  the  four  panes  Instrumented 
wltb  strain  gages  to  indicate  stress  In  the  center  of  the  pane.  The  data 
plotted  In  Pigs.  D.7  and  D.8  give  the  following  relations 

Sfa  =  4.4?  (D.2) 

Sm  =  0.82??  (D.3) 

The  bending  stresses  ere  about  65$  of  those  given  by  Levy  (1942)  while 
the  membrane  stresses  ere  65$  to  80$  of  Levy's  values.  It  Is  Improbable 
that  these  differences  can  be  explained  by  the  fact  the  Poisson's  ratio 
for  the  glass  was  ebout  0.23,  while  Levy  used  v  =  0,316.  Becsuse  of  the 
coefficient  of  4.4  Instead  of  5.91  es  from  tbe  linear  analysis,  stress- 
deflection  results  do  not  coincide  with  the  linear  theory  even  for  small 
deflections . 

Equations  d.2  end  D.3  are  relations  between  the  total  stress  In  the 
center  of  the  plate  and  the  total  deflection  at  that  point.  In  the  ana¬ 
lytical  work  it  will  be  necesssry  to  have  expressions  relating  stress  ano 
deflect  lexis  In  the  flrat  mode  only.  According  to  the  linear  theory  of 
Appendix  A  the  total  static  stress  Is  0.897  of  the  first  mode  stress;  the 
total  static  deflection  Is  0,976  of  the  first  mode  deflection.  We  will 
assume  that  these  factors  from  the  linear  theory  are  velld  for  bending 
but  not  for  membrane  stresses;  that  Is,  the  total  membrane  stress  is 
contributed  by  the  first  mode.  Then 
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NONDIMENSIONAL  BENDING  STRESS 


NONDIMENSIONAL  CENTRAL  DEFLECTION, £ 

T3-«0«5-ZI 

F'G.  D.7  BENDING  STRESS  AS  A  FUNCTION  OF  CENTRAL  DEFLECTION 
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NQNOIMENSIONAL  MEMBRAHE  STRESS  ,  S 


NONDIMENSIONAL  CENTRAL  DEFLECTION ,  ( 

TA-tOM-'fl 

FIG.  D.8  MEMBRANE  STRESS  AS  A  FUNCTION  OF  CENTRAL  DEFLECTION 
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sbi  *  <ri57,<0-878)s>  ■  "■85>  <D-4) 

S  =  0.82  (0.976  ?J2  =  0.7q?,2  (D.S) 

where  the  subscript  1  refers  to  the  first  mode.  In  some  cases  we  will 
want  to  relate  first  mode  stresses  to  total  deflections.  The  appropriate 
relations  are 


Sbl 

-  o!s975  '  <•»?. 

(D.8) 

Sml 

-  0.82?2 

(D.7) 

or  combining, 

Si  = 

=  4,95  (1  +  0.167?) 

(D.8) 

Stress -Loading  Relations 

The  bending 

of  the  loading. 

and  membrane  stresses  can  be  compared  as  direct 

From  the  previous  data  reductions  we  have  that 

functions 

Q  =*  21.7?  +  2.80?s 

(D.l) 

Sb=  4.4? 

(D.6) 

S  =  0.82?z 

m  s 

(D.7) 

Substituting  we 

obta in 

Q  = 

4.93  Sw  +  0.0329  S  * 

D  D 

(D.9) 

and 

Q  » 

24.0  S  ‘/2  +  3.77  S  V2 
m  m 

(d . 10) 

The  data  of  Bowles  Bnd  Sugarraan  (1952)  and  Freynlk  (1963)1  are  compared 
with  these  equations  in  Figs,  d.9  and  D.10.  The  theoretical  results  of 
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NON  DIMENSIONAL  MEMBRANE  STRESS 


BOWLES  AND  SUGARMAN 

O  '/•  "  PLATE 
•  Vie"  SHEET 
A  '/4  "  PLATE 
A  Ve "  PLATE 

FREYNIK 

□  Vo"  WINDOW  GLASS 


1/2 

0-24.0  SmZ+  3.7/  Sm 
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THEORY  OF  LEVY 
(»  ■  0.3(6) 


40  60  80  100 

NONDI  MEN  SIONAL  PRESSURE  ,  0 


TB-SOI5-J4 


FIG.  D.  10  RELATION  BETWEEN  CENTRAL  MEMBRANE  STRESS 
AND  APPLIED  PRESSURE 


Levy  (1942)  sre  slso  included.  Equstion  D.9  sgrees  fairly  well  with 
Levy's  results;  the  dsts  of  Bowles  snd  Sugarman  sppesr  to  fit  Eq.  D.9 
somewhat  better  than  the  curve  from  Levy.  Tbe  membrsne  streaa  results 
ahow  s  marked  difference  between  the  curves  from  Levy  snd  Eq.  D.10.  The 
experimentsl  dsta  sre  divided  between  the  two  curves. 

All  of  the  stress-loading  dats  sbove  is  for  stresses  from  one-quarter 
to  one-half  the  ultimate  stress.  Therefore,  prediction  of  the  stress  con¬ 
dition  st  failure  would  require  some  extrspolation  from  the  svsilsble  dats. 
Because  the  dsta  do  not  ccrrslste  well  with  the  theoreticsl  or  qusal- 
theo:  etical  curves  at  high  stress  levels,  it  does  not  sppear  that  stresses 
st  failure  can  be  predicted  now.  Thia  situation  ia  unfortunste  because 
it  mesns  that  ultimate  strength  dsts  of  gisss  specimens  cannot  be  uaed 
directly  to  predict  fsilure  of  gisss  psnes. 

Statistics  of  Failure 

Strengths  of  gisss  specimens  snd  window  psnes  sppesr  to  vsry  widely 
from  test  to  test;  auch  variations  are  commonly  observed  with  brittle 
materials.  The  experimentsl  results  hsve  usually  been  reported  in  terms 
of  tbe  normal  distribution  although  the  dsts  msy  fit  some  other  distri- 
button  equally  well. 

Table  D.4  lista  the  coefficients  of  varistion  (standard  deviation 
divided  by  the  mesn)  of  the  bursting  pressure  and  the  centrsl  deflection 
for  the  data  of  Bowles  snd  Sugarman  (1952),  The  vsriations  sppear  to 
incresae  with  the  thickness  of  the  specimens,  ss  expected  for  a  brittle 
materisl.  It  should  be  noted  thst  these  variations  pertsin  to  the  pres- 
rure  snd  deflection,  not  to  the  normalized  quantities,  Q  snd  §.  There¬ 
fore,  the  ststed  vsriations  include  the  effects  of  randomness  in  the 
gisss  atrengtb  and  in  the  gisss  dimensions. 

Additionsl  dsts  on  the  statistics  of  the  strength  of  gisss  sre  listed 
in  Table  D.5.  The  value  given  by  Orr  (1957)  is  for  his  estimates  of  the 
stress  at  failure  in  hia  glaaa  panels.  The  other  dats  sre  from  strength 
testa  on  other  types  of  gisss  samples.  Note  that  these  coefficients  of 
variation  pertain  to  stress  at  failure  wheress  the  dsta  of  Tsble  D.4 
refer  to  applied  pressure  snd  ultimate  deflection.  Because  we  do  not 
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have  s  sure  relation  between  atreaa  and  applied  preaaure  at  failure,  we 
cannot  relate  the  coefficienta  of  variations  of  theae  quantities.  Hence, 
a  knowledge  of  the  atstiaticsl  variation  of  breaking  atresa  doea  not  lead 
to  s  knowledge  of  the  statistical  variation  of  bursting  pressure. 


TABLE  D .4 

STATISTICAL  DATA  ON  FAILURE 
(Data  from  Bowles  and  Sugannan) 


Sample 

Description 

No.  of 
Panels 

Bursting  Pressure 
Coeff.  of  Variation 

w 

Central  Deflection 
Coeff,  of  Variation 

(*) 

l/8"  plate 

40 

17  ,2 

8.6 

3/l6"  plate 

30 

17  .9 

9.2 

1/4"  plate 

30 

25.0 

12.1 

3/s"  plats 

30 

23.7 

13.9 

24  oz  sheet 

30 

14  .0 

7,55 

32  oz  sheet 

30 

13  .9 

7.17 

3/l6”  sheet 

30 

26.8 

10.95 

TABLE  D.5 

STATISTICAL  DATA  ON  BREAKING  STRESS 


Source 

Coefficient  of 
Varistion 

w 

Type  of  Glaas 

Frownfelter  (1959) 

22.0 

-  - 

Orr  (1957) 

19.7 

Plate  and  aolex  psnea 

McKinley  (1964) 

25, 0* 

Flate  and  window 

Haaselman  and  Fulrath  (1966) 

12.7 

Sodium  Boroailicate 

Vslue  recommended  by  McKinley  for  uae  in  estimating  safety  factors 


APPENDIX  E 


DYNAMIC  AMPLIFICATION  FACTOR:  NONLINEAR 


When  plates  deflect  beyond  s  small  deflection  range,  their  apparent 
stiffness  Increases  because  of  membrane  action.  This  change  in  stiffness 
increaaea  the  natural  frequency  and  decreases  the  amplitude  of  motion. 

The  purpose  of  thia  appendix  ia  to  determine  the  magnitude  of  these 
changes  for  the  problem  of  window  response  to  booma .  n. e  solution  ia 
developed  first  for  the  motion  in  the  fundamental  mode  only,  and  then 
estimates  are  made  for  the  effects  of  the  higher  modea. 


Firat  Mode  Only 

The  equation  of  motion  in  the  firat  mode  haa  the  following  general 

form 


(E.l) 


where 

5  ia  the  nondimensional  central  deflection,  w0/h 

k  ia  atiffneas  which  ia  s  function  of  deflection 

0^  ia  the  maximum  amplitude  of  the  nondimenaional  forcing  function 

f(t)  ia  the  temporal  variation  of  that  function. 

From  the  solution  of  Eq«  E.l,  we  want  to  develop  curves  of  maximum 
deflection  aa  a  function  of  the  period  ratio,  ujt/2tt,  and  of  the  loading. 
Hence,  we  muat  chooae  to  plot  deflection  aa  e  function  of  loading  for 
specific  valuea  of  ujt/2tt  or  aa  a  function  of  <dt/2ti  for  specific  loadings. 
The  latter  courae  seems  appropriate  because  we  can  expect  the  deflection 
maxima  to  vsry  monotonies lly  with  load.  The  load  values  chosen  are  thoae 
producing  specific  values  of  static  deflection.  The  next  atepa  are  to 
introduce  the  nonlinear  relation  for  k(§)  and  to  replace  by  the  static 
deflection  it  produces. 

The  load-deflection  relation  in  Eq.  C,10  will  be  uaed  although  it 
pertaina  to  total  deflection,  not  luet  firat  mode  deflection.  This  is  a 
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fsirly  small  approximation  because  the  difference  between  static  totsl 
and  first  mode  deflections  is  only  2.5$,  The  relation  is 

Q  =  A05  +  B^3  =  k(S)5  ( E .2 ) 


so  that  the  equation  of  motion  is 


da*  ^  A 
V  +  - 


dt 


u  +  ? 

HI  Ar 


53)  = 


v(t) 


k^s)Qn 

"  ^TTTf(t) 


(E.3) 


Note  thst  k(5  )/m  =  (A0§  +  3)/m§  is  the  square  of  the  natural 

S  8  S  3 

circular  frequency,  and  Q^/kC^)  is  the  static  deflection  £s  under  the 
load  Qp.  Now  simplify  the  equation  to  the  form 


+  a (l  +  e?3)  =  a(?e  +  e§g3)  f(t)  (E.4) 

where 

a  =  A0/m 
e  =  Bj/A0 

Equation  E.4  waa  aolved  numerically  using  an  Adams  predictor-corrector 
method  for  integration.  The  computer  program  uaed  had  been  previously 
developed  for  a  three -degree-cf -freedom  ayatem  and  reported  by  Bycroft 
(1965).  Tha  program  waa  simplified  slightly  and  changed  to  accept  a 
sharply  defined  N-wave  as  the  loading.  The  deflections  and  velocities 
were  printed  out  after  each  step  in  the  integration.  Positive,  negative, 
and  free  vibration  maxima  were  taken  from  these  listings  by  hand.  The 
results  are  plotted  in  Figs.  12-15. 


The  program  is  lieted  st  the  end  of  this  sppendix.  It  is  written 
in  Algol  for  use  on  a  Burroughs  B5500.  The  compilation  time  is  40  seconds, 
and  the  execution  time  per  abscissa  value  averages  14  seconds.  The  input 
data  ia  entered  on  two  types  of  csrds.  An  example  of  input  is  given  in 
the  comment  cards  st  the  beginning  of  the  program.  The  first  csrd  lists 
the  number  of  abscissa  values  for  which  calculations  are  to  be  made.  One 
csrd  of  the  second  type  is  required  tor  each  abscissa.  On  the  second  type 
of  crrd  the  variables  are 
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XI 

ABSCISSA 

FINAL 


EPS 


§  ,  the  stetic  deflection  under  the  maximum  pressure  in 

the  loading  function 

xt/2it 

time  at  which  celculetion  ia  termineted.  It  should  be 
set  et  some  velue  greeter  than  1.0  such  thet  et  ieect 
one  free  vibretion  maximum  ie  reeched. 
e,  the  coefficient  of  the  nonlinear  term. 


In  the  progrem  the  forcing  function  ie  the  etetic  deflection  times 

a  temporal  function.  If  s  response  to  e  particuler  pressure  level,  Q^, 

is  required,  the  appropriate  vf lue  of  g  cen  be  found  from  the  equetion 

6 

%  ■  A°<5.  +  *S„s) 


Alternetively,  the  32nd  cerd  of  the  progrem  could  be  changed  from 


to 


ALFA  (XI  +  EPS  x  XI  *3)  x  OMEGA  *  2; 

ALFA  (Qp/A)  x  OMEGA  *  2; 


where  the  vslues  of  O  end  A  would  have  to  be  supplied. 


The  progrem  output  lists  the  velues  of  ABSCISSA  and  XI  end  then 
prints  four  columns  of  figures,  heeded  TIME,  DEFLECTION,  VELOCITY,  end 
DEFL/XI .  The  velues  printed  under  these  heedings  are  t/T,  g,  Tdg/dt, 
and  |/?s-  The  letter  is  the  dynemic  emplif icetion  factor. 


Higher  ModeB 

A  rigorous  ene lysis  of  the  contribution  of  higher  modes  is  not  made 
tor  the  large  deflection  esse.  However,  some  estimates  ere  required  of 
the  importance  of  these  higher  modes.  From  Appendix  A  the  contribution 
of  these  modes  to  the  lineer  solution  is  known.  What  is  needed  now  is 
e  meens  for  estimating  changes  in  the  amplitudes  end  frequencies  of  the 
higher  modes.  Such  estimates  will  be  made  under  the  essumption  that  the 
motion  in  the  higher  modes  can  be  merely  added  to  that  for  the  fundamental 
mode.  With  this  assumption,  the  problem  cen  be  treeted  es  thet  of  e  plete 
with  stresses  in  Its  plene. 
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The  deflection  of  a  piste  with  stresses  in  its  plane  can  be  derived 
from  section  93  of  Timoshenko  (1959).  The  result  is 


16Si 


Z 

li—  1 


z 

n=l 

odd 


a 


r  /m2+n2.  2 

mn  C(-p — ) 


in (mTTx/a )aln  fnny/a) 

9  f  m  2Nx  +  n*Nv  +  2 


"rr^a  aD 


2amNayij 


(E.5) 


where  N  ,  N  ,  and  N  are  forces  per  unit  width  in  the  midplane  of  the 
x  y  xy 

plate.  They  are  constants  in  Eq.  E.5.  In  our  dynamic  case  Nx,  Ny,  and 

N  all  vary  over  the  surface  of  the  plate,  oscillate  at  the  fundamental 
xy 

frequency,  and  may  be  compreaaive  or  tensile.  However,  Eq.  E.5  gives  an 

indication  of  the  nature  of  the  effect  of  the  membrane  stresses  on  the 

deflections  in  the  higher  modes.  In  the  central  region  of  the  plate, 

the  membrane  atresaee  are  tenaile  and  the  shear  is  zero.  Therefore  the 

assumption  ia  made  that  N  =  N  and  N  =0. 

x  y  xy 


Then 


16tlo  V  r  ain  (mTix/a)  ain  (nny/a) 

™  ctl  n=l  »(»*+"*)/«*  [(m’+n’Va5  +  Nx/tt5d7  (E’6) 


w,n  odd 


The  momenta  are  computed  from 


,,  _  /92w 

M  =  -D  (t— y  +  y  w),  etc. 
x  ' 3x 2  3y  *  * 


(E.7) 


The  momenta  are  then 


00  00  .  , 
16q0a*(l+u)  V1  V"  H(m,n)  ain  (mTTx/a)  ain  (nny/a)  t>y 

“x  =  - V -  i-  L  mAltm^n4)  +  N  aV^oJ -  (E>8) 


m=l  n=m 
odd  only 


where 


H(m,n)  =  1/2  for  m  =  n 

=  1  for  m  =  n 

Evidently  each  term  in  the  series  for  moment  ia  reduced  by  the  factor 


1  +  N^a  a/[n5D(raa+n2) J 
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from  the  value  it  would  have  if  there  were  no  membrane  streasea.  Using 
thia  factor  a  general  expreaaion  can  be  developed  to  relate  reduction  of 
amplitude  in  the  m,nth  mode  to  reduction  in  the  first  mode.  Let  R  (m,n) 
be  the  ratio  of  amplitude  for  the  nonlinear  case  to  that  for  the  linear 
case.  Then 


_ 1 _ _ 

1  +  2[l/"RE(l>l)-Ij/(mz+n2) 


(E.9) 


Aa  an  example  let  ua  aaaume  that  the  amplitude  in  the  fundamental  mode 
waa  reduced  by  30$.  Then  by  how  much  are  the  higher  modes  reduced? 

The  results  for  this  example  are  in  Table  B.l. 


TABLE  E . 1 

AMPLITUDE  REDUCTION  FOR  HIGHER  MODES 


Mode 

N^a  2/rr2D(m2+n*) 

1,1 

0.70  (aaaumed) 

0.428 

1,3  and  3,1 

0.92 

0.086 

3,3 

C.95 

0.048 

3,5  and  5,3 

0.97 

0.025 

The  values  in  the  table  indicate  that  the  amplitudes  of  the  higher 
modes  are  reduced  by  one-fifth  (or  lesa)  as  much  aa  the  first  mode. 

The  natural  frequencies  are  also  altered  by  large  deflectiona.  The 
estimate  of  the  frequency  change  is  made  on  the  same  basin  as  the  arr.pli 
tude  estimate  above.  That  is,  the  modea  are  considered  separately  and 
the  in-plane  stresses  caused  by  the  fundamental  mode  are  the  aource  of 
the  variationa  from  the  linear  theory.  The  approach  used  is  to  calculate 
the  strain  energy  and  kinetic  energy  and  equate  them  to  find  the  natural 
frequencies . 

Again  taking  the  reaulta  of  paragraph  93  of  Timoshenko  (1959), 
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the  deflection  Is  assumed  to  have  the  form 


w  =  Z  ^  <p  aln(mnx/a)  sln(nrry/a)  (E.10) 

m=l  n=l  ”'n 
odd  only 


Thei.  the  strsln  energy  Is 


V 


CO  CO 


E  E 


<p2  [(m2+n2)2  + 

r  m,n  ' 


(maN  +n 2N  +2mnN  ) 
x  y  xy 


(E.H) 


The  kinetic  energy  la 


T 


YhB  2 
8g 


E 


m=  1 


E  (jf  («."))■ 

n=l 


(E.12) 


where 

Y  la  the  unit  weight  of  the  plate 
h  Is  the  piste  thickneas. 

By  equating  the  maximum  values  for  esch  mode  of  T  and  V  and  assuming  that 
the  motion  Is  sinusoidal,  we  can  obtain  the  following  equation  for  fre¬ 
quency: 


m  =  — y  ^J~gD/YL  [(m2+n2) 2  +  (m*N  +n2N  +2mnN  )a  2/tt 2D ] x/ 2 

v  y  n  8  x  y  xy 

Again,  let  N  =  N  ,  N  =0.  Then 
’  x  y>  xy 

«m  n  =  pVgD/YL  ( m 2+n 2 )  [l+a  2Nx/TT2D(m2+n2) ] */2  (F.13) 

The  fsctor  In  Eq.  E.13  rsiaeti  to  the  one-half  power  containa  the  effect 
of  the  ln-plane  atressea.  As  might  be  expected,  this  is  the  same  factor 
which  entered  the  expreaaion  for  the  change  in  amplitude  aa  a  function 
of  in-plane  stresses.  If  the  change  in  frequency  ia  known  for  the  funda¬ 
mental  mode,  then  Eq.  E.13  permits  a  computation  of  the  change  for  all 
modes.  That  is,  if  B(m,n)  ia  the  ratio  of  frequencies  in  the  nonlinear 
and  linear  cases,  thon 
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(E.14) 


BE(m,n)  =*  "Y1  +  2[BE2(l,l)-l]/(ra2+n2) 


Some  sample  results  with  Eq.  E.14  sre  provided  in  Tsble  F.2.  Note  that 
the  aame  valuoa  of  membrsne  stress.  N  ,  are  used  in  Tables  E.l  snd  E.2. 
There  is  evidently  very  little  change  in  the  frequency  of  the  higher  modea. 


TABLE  E.2 

FREQUENCY  INCREASE  FOR  HIGHER  MODES 


Mode 

Frequency  Ratio 

Nxa  2/ir2D(m2+n2) 

1,1 

1.20  (assumed) 

0.428 

1,3  snd  ),1 

1.04 

0,086 

3,3 

1.025 

0.048 

3,5  snd  5,3 

1.01 

0.025 

Tile  two  snalyses  above  hsve  shown  that  the  contribution  of  the  higher 
mrdea  to  deflections  snd  stresses  sre  very  little  affected  by  the  non- 
linesrities.  Therefore  it  is  suggested  that,  for  purposes  of  estimation, 
the  linear  contributions  of  the  higher  modes  be  assumed  to  be  vslid  for 
lsrge  deflections  also. 


PROGRAM  NONLINEAR  WINOOW 


BEGIN  COMMENT  LYNN  SEAMAN  EXT  3567 

PROGRAM  OF  KRIEBEL  WAS  MODIFIED  BY  SEAMAN  (JAN  67)  TO  CALCULATE  THE 
RESPONSE  OF  A  1  DOF  WINOOW  TO  AN  N-WAVE  PRESSURE  PULSE.  PROGRAM  USES 
RUNGEKUTTA  METHOD  AS  A  STARTER  ANO  CONTINUES  WITH  ADAMS  PREDICTOR- 
CORRECTOR  METHOD. 

OATA  IS  ENTERED  ON  TWO  GROUPS  OF  CARDS!  ONE  FOR  ITER.  ONE  FOR  XI. 
ABSCISSA*  FINAL  AND  EPS.  ITER  IS  THE  NUMBER  OF  CASES  TO  BE  HANDLED  ANO 
EQUALS  THE  NUMBER  OF  DATA  CARDS  TO  FOLLOW.  XI  IS  THE  STATIC  OEFLECTION, 
ABSCISSA  =  0MEGATAU/2PI *  FINAL  IS  THE  LAST  TIME  FOR  WHICH  A  CALCULATION 
IS  MADE*  AND  EPS  IS  THE  COEFFICIENT  OF  THE  NONLINEAR  (CUBIC  >  TERMI 
COMMENT  SAMPLE  DATA  CARDS 

1* 

3.0 »  1 .6*  1.5.  0.129.  END  SAMPLE! 

REAL  INTIME. HZERO. INITIAL. FINAL. PRINT. PELB.A-.SB.OELTAT. 

ABSCISSA. XI t OMEGA. EPS .DELTA *  ALFA. DX I ! 

INTEGEK  II. JJ. SIZE. NFIN. ITER. ICI 

REAL  ARRAY  YINITI AL. YF INALC 1 ! 2 3 ) 

FILE  CR  (l.lS)i 

FILE  LP  4(1. 1S>) 

FORMAT  OUT  CLOCK  (X36. "CALCULATION  TIME  OF  PROGRAM  =".F9.2.X3. 

"SECONDS". X36) ) 

INTIME  ♦  TIME  (1)1 

SI2E+2)  NFIN*1 )  DELTAT*1.0)  INITIAL*0.0)  RELB*0.001)  A8SB*0.00I) 

READ  (CR./.ITERH 

FOR  IC*I  STEP  1  UNTIL  ITER 
DO  BEGIN 

0MEGA*0)  ALFA*0)  DELTA*0 )  EPS*0) 

FOR  11*1  STEP  1  UNTIL  2  DO  YINITJALC II 1  ♦  0) 

REAP  (CR. /.XI. ABSCISSA. FINAL. EPS)  I 

H2ER0  ♦  1. 0/(60. OXABSCISSA)!  PRINT  ♦  HZEROI 

OMEGA  *  ABSCISSAX6.2B32) 

ALFA  *  (XI+EPSxXI*3)X0MEGA*2l 

begin 

REAL  TESTIME1.TESTIME2.TPI 

REAL  ARRAY  INFCOSNFIN]) 

REAL  PROCEDURE  FINPUT  (T.FIN)) 

VALUE  T) 

REAL  T  J 

REAL  ARRAY  FINtOll 

BEGIN  COMMENT  SECOND  ORDER  INTERPOLATION  OF  INPUT  ACCELERATION. 

USE  OF  FORWARD  DIFFERENCES) 

REAL  IR.R.Rl »R2l 

INTEGER  II.II1.II2) 

IF  TSTESTIMEI  THEN  BEGIN 
IR  ♦  T/DELTAT) 

II  »  ENTIER  (IR) ) 

W  »  TO.TTt 

IF  ABS ( R ) 58-06 " THEN 
FINPUT  *  FINCII3 

ELSE  BEGIN 

IF  TCTESTIME2  THEN  BEGIN 

III  ♦  11*1) 


102 


1 12  *  III  +  II 

HI  *  I-RI 

R2  *  2-RI 

F INPUT*  O.S*Rlx{R2xFINCII 3-RxFINCII23»+RXR2*FINCIII Dl 

ENO 

ELSE 

IF  KTESTIMEI  THEN  BEGIN 
III  »  II+II 

HI  *  I-RI 

H2  *  2-RI 

FINPUT  *  RlxFINC II 3  ♦  RxFINClIl3l 


ENO 

ELSE 

FINPUT  *  01 
ENO  OF  FINPUT* 

PROCEOURE  FUNCT  (TTEMP, YTEMP.OERIV ) I 

REAL  TTEMP* 

REAL  ARRAY  YTEMP . OERI VC  1 3 1 

bEGlN  COMMENT  COMPUTATION  OF  DERIVATIVES* 

HEAL  Q12  >  Q 1 3 1 

INTEGER  II 

G12  *  YTEMPL  I  .lx  YTEMPC  1  31 

Q13  *  YTEMPC I 3xQ12» 

OEHIVti:  *  YTEMPC 2  31 

OERIVC23  *  -(OELTA  XYTEMPC23  ♦  0MEGA*2  x(YTEMPCl3  ♦  EPS  XQI3>>  ♦  ALFA 
x  FINPUT (TTEMP , INFC  *3) I 
ENO  OF  FUNCT I 

REAL  ARRAY  KFORAOAMSCOJ3.IJ303.  YINCFORAOAMSC I : 303 I  COMMENT  GOES 
BEFORE  AO AMS I 

PROCEOURE  A0AMS1SIZE,  HZERO*  INITIAL,  FINAL*  PRINT.  RELB*  ABSB*  AOAMSOOI 

YINITIAL*  YFINAL.  FUNCT >  I 

VALUE  SIZE.  HZERO.  INITIAL.  FINAL.  PRINT.  RELB.  ABSB  * 

INTEGER  SIZE  I  REAL  HZERO.  INITIAL.  FINAL.  PRINT.  RELB.  ABSB  I 
REAL  ARRAY  YINITIAL.  YFINALC 1 3  I 
PROCEOURE  FUNCT  I 

BEGIN  COMMENT  AOAMS  VERSION  OF  APRIL  I.  1954  I 
OWN  INTEGER  I.  J.  N  I 

OWN  REAL  X.  H.  RELTEST.  ABSTEST*  FACTOR  I 

REAL  H24.  LB.  BOUND.  T.  TEMPX*  TEMPH  I  BOOLEAN  TEST  » 

REAL  ARRAY  Y.  FCO!4.i:SIZE3.  E.  YPUJSIZE3  I 
LABEL  STARTI,  START2,  STARTS'  MARCH.  LASTSTEP.  RETURN  I 
OEFINE  LOOPI  =  FOR  I  *  1  STEP  i  UNTIL  N  00  4  I 

FORMAT  OUT  MSSG  (X4."THE  STEP  SIZE  IS  NOW". EIB. 1 1 > . 

SIMGY 

(X30."**E0UATI0N5  CANNOT  BE  SOLVEO  WITHIN  THE  GIVEN  ERROR  BOUNOS**", X301 


.TITLE 

(XIO. "SOLUTION  OF  NONLINEAR  WINOOW  M0TI0N"//VXI0*"ABSCI5SA  =”,Fs.4, 

X5. "XI  F8.4//XB. "TIME". XS. "OEFLECTION" .X3." VELOCITY", X3."0EFL/XI"//I. 

FRMT  (X5»F8.4.X4,Fa.4,X4,Fft.4,X3,F«.i4»  J 
LIST  LISTI  (ABSCISSA, XI) I 

PROCEOURE  RUNGEKUTTAt YOLO.  FOLO,  YNEW,  FUNCT)  I 

REAL  ARRAY  YOLO,  FOLO.  YNEWC 1 3  I  PROCEOURE  FUNCT  » 

BEGIN 
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OEFINE  K=KFORAOAMStt.  YINC=YINCFQRA0AMSBI 
REAL  INC.  H6  I  INTEGER  L  » 

L  *  0  I  H6  *  H/6.0  ) 

LOOPI  KC0.I3  *  FOLOCI3  I 
FOR  INC  ♦  H/2.0.  INC.  H  00 

BEGIN  LOOPI  YINCCI3  *  YOLDC  I  3+INCxkCL.  I  ]  I 
L  *  L  +1  )  FUNCT (X+INC .  YINC.  KCL»*3) 

ENO  I 

LOOPI  YNEWC  1 1  *  YOLO! I3+H6X(kC0.I 3  '/2.0x(KC  1 . I  3*kC  2. 1  3>  +KC3.n)l 
X  ♦  X  +K 

ENO  RUNGEKUTTA  ) 

BOOLEAN  PROCEDURE  ERRTEST ( YP»  YC.  E?  » 

REAL  ARRAY  YP»  YC .  EC  11  I 
BEGIN  REAL  YCI.  El  I  LABEL  RETURN  I 
ERRTEST  *  FALSE  ) 

LOOPI  BEGIN  ECI]  *EI*  ABS(YPCI3-(YCI*YCCI3))  ) 

IF  El  <  ABS(YCI)XRELTEST  THEN  ECU  ♦  EI/ABS(YCI) 

ELSE  IF  El  <  ABSTEST  THEN  EC  13  ♦  EIxFACTOR 
ELSE  BEGIN  ERRTEST  ♦  TRUE  )  GO  TO  RETURN  END  ) 

END  I 

RETURN:  ENO  ERRTEST  I 
COMMENT  INITIALIZE  ) 

WRITE  (LPCPAGE  3) ) 

WRITE  (LP. TITLE. LISTl) ) 

N  «•  SIZE  )  LOOPI  YC0.I3  ♦  Y INITIAL! I  3  » 

X  *  INITIAL  I  FUNCT  (X.  YC0.*3.  FC0»*3)  ) 

0 XI* Y INITIAL! 1 3/ XI I 

WRITE  (LP.FRMT.INITIAL.LOOPI  YlNITIALC I  3.0X1 ) I 
BOUNO  *  INITIAL+PRINT  I 
IF  CTEST*  ABSBFO)  THEN 

BEGIN  RELTEST  »  14.2xRELB  »  ABSTEST  *  14.2XABSB  » 

FACTOR  *  RELB/ABSB  I  LB  *  RELTEST/200.0  ) 

H  *  PRINT  I  FOR  H*0.5xH  WHILE  H  >HZERO  00  I 

END 

ELSE  BEGIN  rl*HZERO)  IF  H>PRINT  THEN  PRINT*H)  GO  TO  STARTS  END  ) 
COMMENT  RUNGE-KUTTA  STARTING  METHOO  ) 

START1 :  H  ♦  2.0XH  ) 

IF  X+H2FINAL  THEN  BEGIN  J  *  0  )  GO  TO  LASTSTEP  ENO  I 
RUNGEKUTTA <YCO»* 3.  FC0.*3.  YP.  FUNCT)  I  X  ♦  X-H  I 
START2J  H  ♦  O.SxH  ) 
starts:  if  x+h=x  then 

BEGIN  OXI  ♦  YC0.13/XII 

WRITE  (LP.SINGY) I  WRITE  (LP.FRMT.X. LOOPI  YL0.I3.0XI) I 
GO  TO  RETURN 
ENO  I 

RUNGEKUTTA (YC  0**3.  FC0.*3.  YCi**3.  FUNCT))  FUNCT(X»  YCl.*3.  FCl.*3>) 
RUNGEKUTTA (YC 1 . *  3.  FCl,*3»  YC2.*3»  FUNCT)) 

IF  TEST  THEN  IF  ERRYESTUP.  y!2.*3.  E)  THEN 
BEGIN  LOOPI  YPCI3  *  YCI. I  3  ) 

X  ♦  X-H  I  GO  TO  START2  ) 

ENO  ) 

FUNCT(X.  YC2.*3.  FC 2 . *  3 )  ) 

WRITE  (LP.MSSG  .H)i 

IF  X+H2FINAL  THEN  BEGIN  J  ♦  2  »  GO  TO  LASTSTE?  ENO  > 

RUNGEKUTTA (YC 2. *3.  FC2**3»  YC3»*3*  FUNCT) I  H24  *  H/24.0  I 
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COMMENT  CHECK  STARTING  VALUES  FOR  PRINTING  I 
T  *  X-3.0XH  1 
FOR  J  ♦  I r  2*  3  DO 
BEGIN  T  »  T+H  I 
IF  T=BOUND  THEN 

BEGIN  BOUND’BOUND+PRINTI  DXI  *-  YCJfI3/XII 
WRITE  < LP r FRMT > T • LOOPI  YCJfI3fOXI>  END 
ELSE  IF  T>B0UNU-H24  THEN 

BEGIN  TEMPX  v  X  I  TEMPH  *  H  I 
X  ♦  T-H  1  H  *•  BOUND-X  I 

RUNGEKUTTAt YCJ-If*3f  FCJ-1f*3f  YPf  FUNCT)  I 
DXI  ♦  YPCI 3/XI I 

WRITE  (LPfFRKTfXfLOOPI  YPCI3fOXI?) 

X  ♦  TEMPX  I  H  ♦  TEMPH  I  BOUND  *  BOUNO+PRINT 


END  I 

END  I 

COMMENT  ADAMS  MARCHING  METHOD  » 

MARCH:  FUNCTIXr  YC  3*  *  3f  FC3f*3)  I 

IF  X+H2FINAL  THEN  BEGIN  J  f  3  I  GO  TO  LASTSTEP  ENO  I 
LOOPI  YPCI]  *  YL3' I  3  +H24x(55.0*FC  3r  I  3  -5<?.0xFC2H3 
+37.0xFCIfI3  -9.0xFC0fI3)  I 
X  +  X  +H  I  FUNCTUf  YPr  FC4**3>  I 
LOOPI  YC4*I3  ♦  YC3»I3  +H24x(9.0xFC4f I  3  +I9.0xrt3Flj 
-5.0xFC2fI3  +FC I r I  3 >  l 
IF  TEST  THEN  IF  ERRTEST<YP>  YC4f*3f  E>  THEN 

BEGIN  LOOPI  BEGIN  YE0fI3  +  YC3fI3  J  Ft  Or  1 3  *■  FC3»I3  END  I 
X  *  X  -H  I  H  *  0.5XH  I  GO  TO  STARTI  I 
END  ) 

IF  X=BOUND  THEN 

BEGIN  BOUND+BOUNO+PRINT  )  DXI  *  YC4f  I 3/XI I 
WRITE  (LPfFRMTfXfLOOPI  YC4fI3fDXI>  END 
ELSE  IF  X>B0UND-H24  THEN 

BEGIN  TEk.°X  *  X  I  TEMPH  ♦  H  I 


X  *  X-H  )  H  ♦  BOUND-X  I 

RUNGEKIJTTA<YC3f*3f  FE3f*3f  YPf  FUNCT)  I  DXl  *  YPC13/XII 
WRITE  JLPfFRMTfXfLOOPI  YPCI3fDXI)I 
X  *  TEMPX  I  H  *  TEMPH  I  BOUND  *  BOUNO+PRINT 
ENO  ( 

LOOPI  BEGIN  ' YC3* I  3  ♦  YC4fI3  I  FC  0  » I  3  *  Ft  1 f 1  3  I 
Ft  I  f  I  3  *  Ft  2' I  31  Ft2F I  3  *  Ft  3 f I  3 1  FC3fI3  *  FC4fI]  I 
END  » 

comment  can  INTERVAL  H  BE  DOUBLED  ) 

IF  H+H>PRINT  OR  NOT  TEST  THEN  GO  TO  MARCH  t 
LOOPI  BEGIN  IF  EC  1 3  >  LB  THEN  GO  TO  MARCH  END  I 
LOOPI  BEGIN  Yt  0  f I  3  *  YC3fI3  I  Ft  0  f I  3  +  FC3fI3  END  I 
H  +  2.0XH  I  GO  TO  STARTI  I 

LASTSTEP:  H  +  FINAL-X  I  PUNGEKUTTA t YC Jf*3f  FLJf*3f  YFINALf  FUNCT)  I 
DXI  *  YFINALC I 3/XI I 

WRITE  (LPfFRMTfXfLOOPI  YFINALC I 3fDXI ) ) 

RETURN:  end  ADAMS  I 
TESTIMEI  ♦  NFINXDELTATl 


ICC  T  I  Mr 


INF103  ♦  1.0)  INFt I  3  *  -1. 01 

ADAMS  (SIZEfHZEROfINITIALfFINALf PRINT fRELBfABSBfYINITIALf YFINALf FUNCT)) 
ENO  I 
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WRITE  (LPCDBLJH  WRITE  (LPCDBUM  WRITE  (Lp.CLOCK.  (TIME  <l)-INTIME)/60)  I 
END  I 

CLOSE  ICR. RELEASE)* 

END. 

?QATA  CR 

1. 

3.0.  ).8.  1.5.  0.129. 

?ENQ  OF  DECK 
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